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Abstract 

We investigate the thermodynamics of a thermal field theory in presence of 
both a baryon and an isospin chemical potential. For this we consider a 
probe of several D7-branes embedded in the AdS-Schwarzschild black hole 
background. 

We determine the structure of the phase diagram and calculate the rel- 
evant thermodynamical quantities both in the canonical and in the grand 
canonical ensemble. We discuss how accidental symmetries present reflect 
themselves in the phase diagram: In the case of two flavors, for small densi- 
ties, there is a rotational symmetry in the plane spanned by the baryon and 
isospin density which breaks down at large densities, leaving a Z4 symmetry. 

Finally, we calculate vector mode spectral functions and determine their 
dependence on either the baryon or the isospin density. For large densities, 
a new excitation forms with energy below the known supersymmetric spec- 
trum. Increasing the density further, this excitation becomes unstable. We 
speculate that this instability indicates a new phase of condensed mesons. 
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1 Introduction 



Generalizations of the AdS/CFT correspondence [1,2] have proven to be very 
useful in describing strongly coupled QCD-like theories. In particular, the AdS- 
Schwarzschild black hole background provides a gravity dual for four-dimensional 
N = 4 SU(N) SYM theory at finite temperature [3]. This is expected to be 
relevant for describing some aspects of the strongly coupled quark-gluon plasma. 
Hydrodynamic methods based on gauge/gravity duality have been applied very 
successfully to describing transport processes [4-21]. 

A further important ingredient for obtaining gravity duals of field theories 
similar to QCD is the addition of flavor degrees of freedom via the embedding 
of D7-brane probes [22]. A D7-brane probe embedded in the AdS Schwarzschild 
background displays an interesting first-order phase transition between embed- 
dings which do not reach the black hole horizon and those who do [23-26]. The 
first kind of embeddings is often said to be a "Minkowski embedding" while the 
second is called a "black hole embedding". The phase diagram is modified in 
presence of a U(l) baryon chemical potential given by a non-zero vev for the time 
component of the gauge field on the brane [27-30]. The effect of a finite baryon 
chemical potential is also studied in the dS^ background, e.g. [31], and in the 
Sakai-Sugimoto model, e.g. [32,33]. 

The spectral functions for the D7 system, calculated according to the holo- 
graphic prescription of [34,35], also display a rich structure, as investigated for 
instance in [36-42]. In particular, as shown in [43], in presence of a baryon chem- 
ical potential, the vector quasiparticle resonances approach the supersymmetric 
meson mass spectrum of [44] for large values of the quark mass or small values of 
the temperature. At small values of the quark mass or large values of the tem- 
perature, the frequencies of the resonance peaks increase with rising temperature. 
Moreover it was shown in [43] that in presence of an SU(2) isospin chemical po- 
tential given by a non-zero vev for the time component of the gauge field on a 
probe of two coincident D7-branes, the vector spectral function displays a triplet 
splitting of the resonance peaks. 

In the present paper we determine the structure of the phase diagram in pres- 
ence of both a J7(l) baryon and an SU (2) isospin chemical potential, of the form 

:W(o _ o o- <-) 

We thus restrict ourselves to an isospin chemical potential in the c 3 direction in 
flavor space. This restriction was suggested and used also in [45]. The values of 
pP and fi are given in terms of fields which depend on the holographic radial 
direction. Some features of this setup are generalized to a general small number 
Nf <C N c of flavor probe branes. For the discussion of the phase structure we 
restrict to Nf = 2. 

Using a suitable linear combination of the flavor components, the action deter- 
mining the brane embedding splits into a sum of Abelian contributions. Moreover 
the choice of (|1.1|) leads to a number of accidental symmetries which we list in 
detail. We discuss the case of arbitrary Nf as well as Nf = 2. In addition to 
a reflection and a permutation symmetry, there is also an approximate O(Nf) 
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symmetry if the baryonic and isospin densities, parametrized by d and d , are 
small. 

We investigate the phase diagram and calculate thermodynamic quantities 
such as the free energy (grand potential), entropy, energy and speed of sound 
in the canonical (grand canonical) ensemble. Generically, a special case arises 
when the normalized baryon and isospin densities coincide, d 1 = d B . The relation 
Id 1 ] = \d B \ determines fixed lines of the permutation and reflection symmetry 
transformations. 

In the canonical ensemble at vanishing isospin density, there is a phase tran- 
sition between black hole embedding^] below a critical baryon density [27]. At 
both finite baryon and isospin density, the according critical point becomes a crit- 
ical line which displays the approximate 0(2) symmetry up to a small derivation 
which we quantify. We proceed by investigating a new first order phase transition 
in the density plane in the limit of large quark mass over temperature ratio. The 
phase diagram clearly displays the symmetry fixed lines \d J \ = \d B \ at which the 
chemical potentials develop a discontinuous gap, which marks the first order phase 
transition. This is similar to what is observed in two-color QCD [46, fig. 1]. In 
the four different regions of the phase diagram, the larger of the two densities 
measures the degrees of freedom in our setup and contributes linearly to the free 
energy. The smaller density contributes nonlinearly and may be interpreted as a 
electrostatic-like charge induced by the non-Abelian charges of the quarks. 

In the grand canonical ensemble, the sum fi B + /j, 1 is relevant. For fi B + 
fj, 1 < M q there is a black hole/Minkowski phase transition which is absent for 
fi B + [i 1 > M q , with M q the bare quark mass. As we discuss in detail, a new 
phase becomes visible in the grand canonical phase diagram. This is related to 
the fact that it is not possible to find a bijective map between the canonical and 
the grand canonical ensemble in the limit of large mass over temperature ratio and 
of chemical potentials of the order of the quark mass. This implies an instability, 
e. g. [47, chap. 3]. 

This motivates us to calculate the spectral function for the gauge field fluctu- 
ations in the canonical ensemble at either finite baryon or finite isospin particle 
density. We identify resonance peaks of the spectral functions with mesonic quasi- 
particles. For finite baryon density, it was shown in [43] that for ji B ~ M c , with 
M c the renormalized quark mass in the plasma, the spectrum of these particles ap- 
proaches the super symmetric spectrum M n /(2irT) = m\J{n J r l)(n + 2)/2, with 
m = 2M q /(\/ r \T) found in [44]. The quasiparticle resonances become very narrow 
in this case. Here, however, we investigate the case of large baryon density and 
find that these resonances become wide, i.e. unstable, in this limit. 

In the case of finite isospin density, we additionally find a new excitation. This 
new resonance develops at an energy below M = 2nTm. The new resonance be- 
comes the lowest vector mesonic excitation of the system. It becomes unstable at 
a critical isospin chemical potential. We interpret this as vector meson condensa- 
tion, similar to [48]. Moreover we discuss the structure of the phase diagram in 
(fj, 1 , Mq/T) and identify the new phase. We relate the appearance of this phase 

1 Note that as discussed in [28], this is really a transition between a black hole and a mixed 
phase, since the lower energy black hole phase has been shown to be unstable. 
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to the instability discussed above. 

We also calculate the effective baryon diffusion coefficient D within hydrody- 
namics from the membrane paradigm. We find that D displays the hydrodynamic 
crossover reminiscent of the phase transition present at small densities. The posi- 
tion of the crossover displays a Z4 symmetry in the (d B , )-plane. 

The paper is organized as follows. In section [2] we introduce the setup of probe 
D7-branes in presence of both a baryon and an isospin chemical potential, paying 
particular attention to the accidental symmetries. We obtain numerical results 
for the chemical potentials and the condensate as functions of the quark mass. 
In section [3] we study the phase diagram and thermodynamic quantities such as 
free energy, entropy and energy in the canonical ensemble. In section H] we per- 
form a similar analysis in the grand canonical ensemble, revealing the possibility 
of an instability. The hydrodynamic diffusion constant is studied in section [5l 
In section [6] we study the spectral functions of mesonic quasiparticle resonances. 
For large isospin chemical potentials, we find an instability and identify the cor- 
responding region in the phase diagram. 



2 Holographic setup 

2.1 Background and brane configuration 

We consider asymptotically AdS§ x S 5 space-time. The AdS§ x S 5 geometry is 
holographically dual to the N = 4 super- Yang-Mills theory with gauge group 
SU(N C ). The dual description of a finite temperature field theory is an AdS black 
hole. We use the coordinates of [27] to write the AdS black hole background in 
Minkowski signature as 



with dr2§ the metric of the unit 5-sphere and 

/(e) = 1-4. /(e) = 1 + 4, (2.2) 

where R is the AdS radius, with 

R 4 = Airg s N c a' 2 = 2A a' 2 . (2.3) 

The temperature of the black hole given by (|2.1jl may be determined by demanding 
regularity of the Euclidean section. It is given by 

r = S- < 2 - 4) 

In the following we may use the dimensionless coordinate p = q/qh, which covers 
the range from the event horizon at p = 1 to the boundary of the AdS space at 
p — > 00. 

To include fundamental matter, we embed Nf coinciding D7-branes into the 
ten-dimensional space-time. These D7-branes host flavor gauge fields with 
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gauge group U(Nf). To write down the DBI action for the D7-branes, we introduce 
spherical coordinates {r, Q3} in the 4567-directions and polar coordinates {L,(f>} 
in the 89-directions [27]. The angle between these two spaces is denoted by 9 
(0 < 9 < vr/2). The six-dimensional space in the 456789-directions is given by 

dp 2 + g 2 dQ 2 5 = dr 2 + r 2 dfi| + dL 2 + L 2 d^ 2 

= d£ 2 + g 2 (d9 2 + cos 2 9d4> 2 + sin 2 9dn 2 ) , ^ 

where r = £sin#, g 2 = r 2 + 1? and L = gcos9. Due to the symmetry, the 
embedding of the D7-branes only depends on the radial coordinate p. Defining 
X = cos 9, we parametrize the embedding by x = x{p) an d choose (j> = using 
the 0(2) symmetry in the 89-direction. The induced metric G on the D7-brane 
probes is then 

ds 2 {G) = JL (-£dt 2 + /df 2 ) + ^ ~ X \ + _ Q ^ deX)2 dg 2 + R 2 (l - X 2 )dn 2 . 

[ ' (2 - 6) 

The square root of the determinant of G is given by 



\—G = ^fg 6 ff(l -xyi-x 2 + 6 2 {d Q X? , (2-7) 

where /13 is the determinant of the 3-sphere metric. The embedding function x 
will be determined numerically in section 12.3.31 It depends on the baryon and 
isospin chemical potentials we introduce next. 



2.2 Introducing baryon and isospin chemical potentials 

We introduce baryon and isospin chemical potential as a non-vanishing time com- 
ponent of the non-Abelian background gauge field living on the D7-branes [27], 

H = lim A (p) . (2.8) 

First we consider the case of two different flavors, Nf = 2. The generators of the 
gauge group U (Nf = 2) can be chosen as the three Pauli matrices a a completed 
with the identity a , 

■"-(J!). ^(io). •--(!?). ^G-O- <-) 

For simplicity, we only consider the diagonal representations of the gauge group. 
This is equivalent to rotating the flavor coordinates until our chemical potential 
lies in the third isospin direction. Thus, the field strength tensor F N f = dA N f + 
A f A A f on the D7-branes may be written as 

f%/ =2 = F*y + . (2.10) 

The part of the field strength F B a° corresponding to the U(l) gauge group charges 
the branes equally and therefore induces a baryon chemical potential [27]. How- 
ever, the part F 1 ^ corresponding to the SU(2) gauge group charges the branes 
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differently inducing a isospin chemical potential [49, 50]. The definition above 
may be generalized to arbitrary Nf > 2. Recall that there are (Nf — 1) diagonal 
generators of SU(Nf) which form the Cartan subalgebra. Inspired by the inter- 
pretation that a diagonal generator of SU(Nf) should charge one brane differently 
with respect to all others, we write the diagonal generators as 

i-th position 

A i = diag(l,...,-(iV/-l),...,l) i = 2,...,N f . (2.11) 

For this choice of matrices the first flavor component is treated as the reference 
quantity. For arbitrary Nf the generator of the baryonic part of the gauge group 
U(l) is called A 1 . Thus, we can generalize (12. 10f> to 

1=1 

where again F B = F Jl induces the baryon and F li for i > 2 the isospin chemical 
potential for the i-th flavor component. In our setup, the only non-vanishing 
components of the background field strength are F^q = —Fq^ = 8 p Aq since Aq 
depends on p only. 



Nf 

f?J = f^ + ^f^ = 

i=2 



2.3 The DBI action and equations of motion 

The DBI action determines the shape of the brane embeddings, i.e. scalar fields <£, 
as well as the configuration of gauge fields A on these branes. We consider the 
case of Nf = 2 coincident D7-branes for which the non-Abelian DBI action reads 

Sdbi = -T D7 Str J d 8 ^ Vdet Q^/det (P ab [£~ + E^Q- 1 - S)^E~} + 2^a'F ab ) 

(2.13) 

with 

Q { j = Fj + \2Tra'{<$>\ <5> k ]E kj , (2.14) 

where for a Dp-brane in d dimensions we have fi, v = 0, . . . , 9, a, b = 0, . . . ,p, 
h J = (P +!))••• >d, Ep, v = g pu + B^ u . In our case we set p = 7, d = 10, 
B = 0. The action (|2.13p can be simplified significantly by using the spatial and 
gauge symmetries present in our setup. First we make use of the spatial rotation 
symmetry in our 8,9-directions in order to rotate to the frame in which <£ 9 = 0. In 
this particular frame, all the commutators of our two scalar fields $ 8 , $ 9 vanish 
and thus = <5!-. However, the non-Abelian structure of embeddings and gauge 
fields is still manifest in the pullback appearing in the action (|2.13|) as 

Pab[g^] =9ab + (2nafg l3 {d a &d b & + id a &[A b , &} 

+i[A a ,&]d b & - [A a ,&][A b ,&}) . { " U] 

The only terms coupling the embeddings $ to the gauge fields A are given by 
the commutator terms in equation (|2.15|) . These commutators vanish due to 
the following argument. We are free to define the r° flavor direction to be the 
direction parallel to the non-vanishing embedding $ 8 using the U(l) C U(Nf) 
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gauge symmetry. Now we choose our baryonic gauge field to have only a non- 
vanishing 74Qdx°T component. In the dual field theory, this particular choice 
corresponds to the case in which the baryon charge representation and the mass 
representation are simultaneously diagonalizable. A different choice is possible 
and will most likely change the field theory phenomenology. This applies also 
to the case of non-coincident Dp-branes. Finally we use the remaining gauge 
symmetry SU(Nj) C U(Nf) in order to restrict to an isospin gauge field A along 
the third flavor direction only, i.e. AqcIx^t 3 , without loss of generality. Thus the 
only representations appearing in the background are the diagonal elements t° 
and r 3 of U(Nf = 2), which constitute the Cartan subalgebra. By definition 
these representation matrices commute with each other and thus all commutators 
[$, A] vanish. Therefore all terms coupling the scalars to the gauge fields A 
vanish, such that their equations of motion decouple as well. In this way we can 
consistently truncate the non-Abelian DBI action to 



S DB i = -T D7 J d 8 £ Str ^ | det (Ga° + 2vra , F Ar /= 2 ) |^ 



T D7 Jd s C v^GStr No-° + (2vra / ) 2 G 00 G 44 (f, 



N f =2 
40 



(2.16) 



where in the second line the determinant is calculated. Next we determine the 
square of the non- vanishing components of the field strength tensor, which is by 
construction diagonal in the flavor space, 



N t =2 



^40 



2 



(Fi) 2 + (Fl f 



a + 2FgFi a 3 



(2.17) 

= diag { {Fg + Fg) 2 , (F& - Fg) 2 } , 

where we used Defining the new fields! 

X l = A$ + A I and X 2 = A$-A I , (2.18) 
the square of the field strength can be written ae[l 

(F^ =2 ) 2 = diag {{dpX^ 2 , (d p X 2 ) 2 } . (2.19) 

Notice that the field Xi is the gauge field living on the i-th brane. Inserting the 
metric components, the action (I2.16P becomes 

Sbbi = ~T D7 [d s Z - X 2 ) 




1 - X 2 + Q 2 {d e x) 2 ~ 2(2W)2-i(l - x 2 )(d,X^ (2 . 20) 
+ Jl-X 2 + S 2 (d gX ) 2 ~ 2(2vra') 2 -i(l - X 2 ){d Q X 2 )A . 



2 To distinguish between the physical basis and the basis, in which the fields decouple, we give 
the quantities in the physical basis as A 3 ^ 1 an upper index and the quantities in the decoupled 
basis as Xi a lower index. 

3 Notice that only the time component of the gauge field A B ^ X does not vanish and depends 
only on the radial coordinate p: F 3 Q = d p A B 
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The transformation to the fields X{ decouple the two branes and we obtain a 
sum of two Abelian DBI actions. The Abelian actions are known from the pure 
baryonic case [27] and we therefore can use the ideas given in the pure baryonic 
case to study our setup in which baryon and isospin charges are switched on 
simultaneously. 

Now we generalize this result to arbitrary Nf. Also for this case, we can 
decouple the non- Abelian gauge field into Abelian gauge fields on each brane. For 
simplicity we suppress the according index Nf on the field strength tensor. We 
again start from the DBI action for the D7-branes, 

Sdbi = ~T D7 Jd s Z Str ( A /|det(GA 1 + 2vra / F)|) . (2.21) 

Since this action is diagonal in flavor space, we are able to evaluate the square 
root and the trace directly (for more details see appendix [A]). After a redefinition 
of the fields 

X i = J2 A Q> X i = J2 A Q ~( N f-V A Q ' i = 2,..-,N f , (2.22) 

3=1 3=1 

where Xi is the i-th flavor component of the gauge field A f, the non- Abelian 
DBI action again becomes a sum of Nf Abelian DBI actions, 



Nf 

< 

1=1 

(2.23) 



E \l 1 - X 2 + Q 2 (d eX ) 2 ~ 2(2vra') 2 -^(l - X 2 )(d e Xrf . 



Thus we may consider each brane separately with the gauge fields X{ treated as a 
U(l) gauge field living on the i-th brane. The action for each brane is the same as 
the action for the pure baryonic case [27]. In the following we apply the approach 
of [27]. First we calculate the electric displacement di on the i-th brane, which is 
a constant of motion and proportional to the flavor charge density, 



d i = t/a v ^ =( 27rQ ) T D7V h 3 Q 777 



8{d e Xi) v ' v °" 2/ 

X 



(1 - X *)%X t ( 2 - 24 ) 



1 - X 2 + Q 2 (d e X? ~ 2(2vra') 2 f (1 - X 2 )(d g X^ 

From the relations of the gauge fields (I2.22p . we can read off the relations between 
the conjugate charge densities 

Nf 

d B = d h = E d J > dh = J2 d J~ ( N f ~ l ) di i = 2 T--> N fi ( 2 - 25 ) 
j=l 3=1 
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which for Nf = 2 becomes 



<F =di + d 2 and d 1 = di - d 2 . (2.26) 

We now construct the Legendre transformation of the action (j2.23j) to eliminate 
the fields in favor of the constants di , 



Sdbi = Sdbi — / d 8 £ } X, 



N f 

^ v <5Sdbi 



-T D7 f d 8 e ^V//(l - X 2 )^Jl-X 2 + e 2 (d e x) 2 (2-27) 



8d? 



x 
i=l 



(27ra') 2 r| )7 ^/ 3 (l-x^ 3 



The gauge fields can be calculated from the Legendre transformed action by 
d e Xi = -SS-DBi/Sdi- 

It is convenient to introduce the dimensionless quantity 



di 

2ira'g 3 H T D7 



di = n .\ „ ■ (2-28) 



2.3.1 Accidental symmetries 

The DBI action (|2.27|) has the following discrete symmetries: 

Permutation: di <-» dj for i,j = 1, . . . , iVj : Since cij is the electric displacement 
of the i-th brane and the coincident branes have the same embedding func- 
tion by definition, the substitution di <-> dj interchanges the two branes. 

Reflection: di «-> —di for each i = 1, . . . , iVj : We consider an U(l) gauge field 
with only the time component non-vanishing. Since it does not depend on 
the time, our setup is like an electrostatic situation. To change the sign of 
the electric displacement d{, we must change the sign of the gauge field Xi 
(see equation ()2.24l) ). The sign change in the gauge field may be induced by a 
sign change of the electric charge generating the gauge field X%, which may 
be interpreted as an interchange between particle and anti-particle. The 
symmetry di *-* —di is therefore induced by the invariance of our system 
interchanging particles with anti-particles on the i-th brane. 

If we expand the action (|2,27[) for small di, 

Sbbi = ~T D7 J d 8 e - X 2 )^l-X 2 + P 2 {d P x) 




(2.29) 



(27ra') 2 T I)7/0 6/3 ( i_ x 2 ) 3 



+ ••• 



4 The definition of d in the pure baryonic case [27] differs from our convention by a factor Nf 
such that the dimensionful physical parameters d 1 ' have the same normalization. 
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we observe an approximate O(Nf) symmetry 



(d 1 ,...,d Nf )^(D(d 1 ,...,d Nf ), OeO(N f ). (2.30) 

The approximate O(Nf) symmetry, which consists of the SO(Nf) group and re- 
flections, is the continuous symmetry generated by the discrete symmetries since 
the SO(Nf) generators are real skew-symmetric matrices which can be represented 
by reflection and permutation matrices. 

In the physical basis the effect of the symmetry transformation above is given 

by 



di «-> dj\ 



for i, j / 1 : d h <-> d Ij 

for i = 1, j / 1 : d B <-> d B and d 1 ' <-» - ^ d 1 ' 



N 



f 

'k 



(2.31) 



k=2 



di <-> — 



d h <-> -d h +2^dj and d 1 ^ <-> d 1 ' -2d; j + i. (2.32) 



For Nf = 2 we get 



d\ <-> C?2 : 


d B 


^d B 


and 


r/ 7 


i— > — d 


d\ <-> — d\ : 


d B 


' ^ -d 7 


and 


r/ 7 


i— * —d 


d 2 <-> —di : 


d B 


' ^ d 7 


and 


r/ 7 


^d B 



(2.33) 



These symmetries are also present in 2-color QCD theories with degenerated quark 
mass [46]. Since the transformation matrix to the physical basis is proportional 
to an 0(2) matrix, there is also an induced approximate 0(2) symmetry in the 
physical basis. The results for the case Nf = 3 are shown in the appendix [Bl 

2.3.2 Gauge fields 

We now determine the gauge fields which will allow us to compute the chemical 
potentials. Therefore, we write down the asymptotic form of the gauge field close 
to the boundary, 



A o=^~^f + --- fOT * = (2-34) 



where the coefficients /i 7i and d li are related to the baryon/isospin chemical po- 
tential and the baryon/isospin density n 7i , respectively, 

5 j 

A - q (2.35) 



NfN c V\T 3 ' 
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Obviously, the chemical potential p Ii depends on the corresponding density d 1% . 
To determine, this dependence it is convenient to use the dimensionless quantities 

Xi = Xi, m = — m = \-—. (2.36) 

QH QH VAT 

As described above, the transformed action (|2.27)) may be used to calculate the 
gauge fields and therefore the chemical potentials (see [27] for more details) 



j* = Xi(p = 00) = 2di I dp , fJE2M2M£ (2.37) 



V/(l-X 2 )[p 6 / 3 (l-X 2 ) 3 + 8^] 



The connection to the physical basis is given by 

p 1 = J2^, Ai = X)A /j -(^/-l)^ i = 2,...,N f , (2.38) 

i=l 

which for iVj = 2 becomes 

/2 B = -(fix + /2 2 ) and fi 1 = -{fix ~ h) ■ ( 2 -39) 

Notice that in the Abelian case [27] it is found that in the limit m — > 00 the 
chemical potentials /Xj approach M g for positive dj. This allows us to calculate the 
physical chemical potentials in the limit m — > 00, 

H B ^ M q and ^ i = 2, . . . , N f , (2.40) 

if dj > for j = 1, . . . , Nf. In the case dj < 0, using the accidental symmetry of 
the action di <-» — presented in section 12.3. H we get also a sign change in the 
chemical potential fii <-> — //j. Therefore, if we change the sign of just one di, there 
is at least one isospin chemical potential jl 1 * which does not vanish any more. For 
Nf = 2 we explicitly get 

ft B -> M q and p} -> for di, J 2 > i. e. < d 1 < d B 



> and / u / -> M g for - d x < d 2 < i. e. < d B < d 1 . 

(2-41) 

If we set one density to zero in the case Nf = 2, e. g. d 2 = 0, only the first brane is 
charged. Therefore, the chemical potential p\ approaches M q and /U 2 approaches 
zero in the limit m — >• 00. Using equation (|2.39p . this implies 

fi B -» and // -> for di > 0, d 2 = i. e. < d B = d 1 . 

(2.42) 

This discontinuous step suggests that there is a phase transition between the 
regions < d 1 < d B and < d B < d 1 . We also expect a similar phase transition 
for arbitrary Nf. In section [3731 we discuss this phase transition further. 
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2.3.3 Embeddings 



The equation of motion for the embedding function x can be derived from (|2.27|) 

as 



P 5 ff(i-x 2 ) 
p 3 ffx 



N 



8d? 



Vl-X 2 + P 2 (d pX ) 2 Uy P 6 / 3 (1"X 2 ) 3 



Vi-x 2 + P 2 (d pX y 2 

24 



[Z{l-X 2 )+2 P 2 {d pX ) 2 ]Y J \ 1 + 



p 6/3 (1 _ x 2)3 



i-x 2 + P 2 (d pX f)Y: 



1=1 

d 2 



p 6 P(i-x 2 ) 



<=i W1 + 



8d? 



P 6 / 3 (l-X 2 ) 3 



(2.43) 



In addition to the case without isospin density [27], sums over the Nf different 
densities appear in the equation of motion for the embedding function. Due to 
these sums, the equation of motion cannot be written down in a more compact 
form, as done in the Abelian case [27]. 

To relate the gravity field x to the dual gauge field parameters we consider the 
asymptotic form for the embedding function % close to the boundary, 



??? 



X 



P P° 



(2.44) 



Then the AdS / CFT dictionary relates the coefficients m and c to the bare quark 
mass M q and the quark condensate ($1/)), respectively, 



m 



2M q 



8(#) 



V\NfN c T 3 



(2.45) 



The initial conditions to solve equation of motion (|2.43p numerically are x{p = 
1) = Xo and d p x(p = 1) = 0. We determine the parameter m and c, which depend 
on xoj by fitting the numerical solution to the asymptotic form (I2.44|) . 

Examples for embeddings are shown in fig. [TJ In the case of vanishing particle 
density we observe both black hole embeddings, ending on the black hole horizon, 
and Minkowski embeddings [23]. The latter embeddings do not touch the horizon. 
For any finite value of the particle density there are only black hole embeddings 
[27]. 



2.3.4 Strings from branes 

In the case of pure finite baryon density [27] it was shown that D7-branes repre- 
senting heavy quarks develop a spike near the horizon, since they reach the horizon 
albeit their large value of L at large r. This spike may be interpreted as a bundle 
of strings stretching between the D7-branes and the black hole. This bundle of 
strings is necessary in the presence of a non-vanishing gauge field on the brane. 
Since in our case we have a similar setup, we expect to have this spike, too. To 
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Figure 1: Embedding function L(r) of D7-branes in the AdS black hole background, 
with the dimensionless coordinates L = pcosO — p\ and r = psinf? = pyl — x 2 ■ 



formalize this intuition, we analyze the near horizon limit of the brane embedding 
in more detail as in [27]. We start by writing the Legendre-transformed action as 



(2.46) 

Notice that x = cos 6, which becomes x ~ 1 if the embedding is very near to the 
axis. Therefore, the second term in the square roots can be neglected and we get 



V 3 vol(S 3 ) / ' \ /■ / / , e 2 {d eX ) 



5dbi = — \ , > di /dtde-7=Wl + 





X 2 



dtdg J-gtt(g ee + gee(d e 9?) , (2.47) 



where V3 is the Minkowski space volume and vol(»S 3 ) the volume of the 3-sphere. 
Recognize the fact that the result above can be written as the Nambu-Goto action 
for a bundle of strings stretching in g direction but free bending in the 6 direction, 

Sdbi = -^ 3 vol(,S 3 ) Sn g ■ ( 2 - 48 ) 

Therefore just one of the densities di must be non-zero to make the branes develop 
a spike close to the axis. This is the expected result. According to the discussion 
below equation (|2.13p we consider a consistent case in which the Nf branes are 
coincident. Therefore they have the same embedding function, such that the 
complete stack of Nf D7-branes develops the spike discussed above. 
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Figure 2: The dimcnsionless mass parameter m as defined in equation (|2.45[) versus the 
horizon value xo = hnip^i X 01 the embedding at baryon density d B = 0.5 for the case 
Nf — 2. The five different curves correspond to isospin density d 1 = (black), d 1 — jd B 
(green), d 1 = ±d B (blue), d 1 = \d B (red) and d 1 = d B (orange). 




(a) (b) (c) 

Figure 3: The dimcnsionless chiral condensate c versus the mass parameter m as defined 
in equation ((2~l5)) at baryon density d B = 5 • 1(T 5 (a), d B = 0.5 (b) and d B = 20 (c) for 
the case Nf = 2. The five different curves in each figure correspond to isospin density 
d 1 = (black), d 1 = \d B (green), d 1 = \d B (blue), d 1 = |d s (red) and d 1 = d B (orange). 
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0.0 L 0.0 



(c) (d) 

Figure 4: The baryon (left) and isospin (right) chemical potential divided by the bare 
quark mass M q versus the mass parameter m as defined in equation (|2.45[) at baron density 
d B = 5 • 10~ 5 in the upper figures and d B = 0.5 in the lower ones for the case Nf = 2. The 
five different curves in each figure correspond to isospin density d 1 = (black), d 1 = jd B 
(green), d 1 = \d B (blue), d 1 = \d B (red) and d 1 = d B (orange). 



2.4 Numerical results for the background fields 
at constant baryon and isospin density 

In this section we present numerical results for the background fields at constant 
densities, which will later be used to determine the thermodynamics in the canoni- 
cal ensemble. We first analyze the embedding of the D7-branes. Fig. [2] and [3] show 
the dependence of the mass parameter m on the horizon value of the embedding 
function xo = lim p ^i \ and of the chiral condensate c on the mass parameter m 
for the case Nf = 2 at different baryon and isospin densities. 

Fig. [2] shows the same divergent behavior of the mass parameter as found 
in [27] for the pure baryonic case. Fig. (a) shows that the low temperature 
black hole embeddings for small but non-zero densities mimic the behavior of 
Minkowski embeddings which describe the brane at zero densities. For small 
densities there is a first order phase transition between two black hole embeddings, 
which replaces the first order phase transition between black hole and Minkowski 
embeddings at zero densities [23,25]. Later in section l3~2l we discuss that the 
phase transition between the two black hole embeddings must be replaced by a 
phase transition between black hole embeddings and an inhomogeneous mixture 
of black hole and Minkowski embeddings since the black hole embeddings alone 
are not the stable ground state of this theory. At a baryon density larger than a 
critical density, which value depends on the isospin density, the phase transition 
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disappears and a local maximum appears in the chiral condensate (see fig. [3] (b)). 
This maximum disappears if we increase the baryon or the isospin density and the 
chiral condensate c monotonically decreases as we increase the mass m (see figE] 
(c))- 

Fig. also shows that for finite densities, the chiral condensate c approaches 
a fixed non-zero value in the limit m — ► do, i. e. T — > 0. In fig. [3] (a) this asymp- 
totic value is small due to the small baryon density but also non-zero as in fig. [3] 

(b) and (c). Since the non-zero chiral condensate breaks conformal invariance 
and supersymmetry, these symmetries are always broken in our setup even in 
the limit T — > 0, where the AdS black hole background becomes conformal and 
supersymmetric. This symmetry breaking is induced by the finite densities and 
chemical potentials. The finite densities, which introduce a new scale, break con- 
formal invariance. The non-zero vev of just the time component of the gauge field, 
which describes the chemical potential, breaks Lorentz invariance and therefore 
supersymmetry. 

Next we present the chemical potentials \x B and /r in fig. 0] for the case Nf = 2 
at different baryon and isospin densities. In the limit m — > do, the baryon chemical 
potential \x B approaches M q and the isospin chemical potential [i 1 approaches 
(see fig. [U black, green, blue and red line) for d B > d 1 > as shown in equation 
(I2.4ip . For d B = d 1 (orange line) the baryon and isospin chemical potential both 
approach M q /2 in the large mass limit as shown in equation (I2.42|) . Moreover, fig.[U 

(c) and (d) show that the baryon chemical potential is independent of the isospin 
density except for the discontinuous step at d B = d 1 and the isospin chemical 
potential obviously depends on the isospin density. Fig. U] (a) and (b) show that 
the phase transition is also visible in the chemical potential. For a baryon density 
smaller than the critical density, the dimensionless quantities fx B ^ /M q decrease 
as we decrease the mass m and stays constant at a non- vanishing but small value 
after the phase transition. For a baryon density larger than the critical density, 
the dimensionless quantities [i®! 1 /M q diverge in the small mass limit m — > 0. 

2.5 Numerical results for the background fields 

at constant baryon and isospin chemical potential 

In the calculations above we treated the densities d B and d 1 ' as a independent 
variables since they are the constants of motion. Thermodynamically this means 
that we consider the canonical ensemble. Since we also would like to study the 
grand canonical ensemble with the chemical potentials jl B and as independent 
variables, we must analyze the background fields for constant chemical potentials 
fi B and jl To determine the densities di for a given chemical potential fxi we 
must invert equation (I2.37p . We use a shooting method to obtain these functions 
di(fjLi,m)E 

So far we considered only the case with non-zero densities. We now investigate 
the case with zero densities but still non-zero chemical potentials. As shown 
in [51], in this case Minkowski embeddings must also be included. Since the action 

5 For the exit condition of the shooting method we use that the deviation from a given chemical 
potential flf to a chemical potential p,i(di) calculated numerically by (|2.37p for a chosen di is at 
most 0.5%. 
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only depends on the derivative of the gauge field, a Minkowski embedding with 
constant gauge field solves the equation of motion. Therefore, on these Minkowski 
embeddings the chemical potential can take any value. The equation of motion 
for the Minkowski embedding L is given by [25] 

a \ 3 (\ - 1 ^ drL 
r [ V {t 2 + l 2 y) yiTTW 

where L and r are dimensionless and given by p 2 = r 2 -\-L 2 , r = p sin 9 = pyl — x 2 
and L = pcosO = px- From now on we only use these dimensionless coordinates 
L and r. The asymptotic solution for the embedding function L close to the 
boundary is, 

£ = m + 4 + ---- ( 2 - 5 °) 

The parameters m and c are related to the quark mass M q and the chiral conden- 
sate {ipip) given by equation (|2.45p . Using the initial conditions L(r = 0) = Lq 
and d r L(r = 0) = 0, we may solve the equation of motion (j2.49j) numerically. The 
parameters m and c may be determined by fitting the numerical solutions to the 
asymptotic form (|2.50p . 

In fig. U and [8] we show numerical results for the background fields at 
constant chemical potential. Fig.[5]shows the mass parameter m versus the horizon 
value of the black hole embeddings xo between and 1 and the asymptotic value 
of the Minkowski embeddings Lq between 1 and 2 at different baryon and isospin 
chemical potentials. In these figures we see that there is a region where the 
Minkowski and the black hole embeddings generate the same mass parameter m. 
We will show in section H~2l that in this overlap there is a first order phase transition 
as in the case of zero chemical potential [23,25]. For small chemical potentials 
(see fig. GO (a)), the mass parameter decreases in the black hole phase for xo — ► 1 
as known from [23,25]. However, for larger chemical potentials (see fig. [5] (b)), 
the mass parameter m increases monotonically as we increase xo- I n both cases, 
fig. [5] (a) and (b) demonstrate that the mass parameter depends linearly on the 
the asymptotic values xo an d Lq in a large region. The only non-linear behavior is 
in the region where both embeddings can be constructed. For chemical potentials 
(p B + p I )/M q > 1 (fig. 0(c)), the mass parameter diverges as xo - *• 1 as in the case 
of constant densities. This corresponds to the fact that in the canonical ensemble, 
there are black hole embeddings for all mass parameters m. In the grand canonical 
ensemble, however, there are also Minkowski embeddings. 

Fig. [6] shows the chiral condensate c versus the mass parameter m at different 
baryon and isospin chemical potentials. For small chemical potentials (fig. [6] (a) 
and (b)), we see a similar behavior as in the case of zero chemical potentials [23,25]. 
By increasing the chemical potentials, the spiral behavior of the chiral condensate 
vanishes since the mass parameter increases monotonically as we increase xo- 
However, there is still a region where the chiral condensate is multivalued since 
there are black hole and Minkowski embeddings which generate the same mass m. 
For large chemical potentials (p B + p r )/M q > 1 (fig. El (c)), the chiral condensate 
diverges as we increase m. Since for smaller chemical potentials (p B + p I )/M q < 1 
the chiral condensate in the black hole phase approaches the value for Minkowski 



r A L 



(r 2 + L 2 ) 1 



Vl + (d r L) 2 = 0, (2.49) 
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(a) (b) (c) 

Figure 5: The dimcnsionless mass parameter m as defined in equation (|2.45[) versus 
the horizon value of the black hole embeddings xo — nm p^i X between to 1 and the 
asymptotic value of the Minkowski embeddings Lq = lim r _>o L between 1 and 2 at baryon 
chemical potential (J, B /M q = 0.01 (a), fi B /M q = 0.1 (b) and n B /M q = 0.8 (c) for the 
case Nf = 2. The dotted purple curve corresponds to Minkowski embeddings and the 
five other curves to black hole embeddings with isospin chemical potential [i 1 = (black) , 
// = jfi B (green), // = \ii B (blue), [i 1 = j[i B (red) and // = \i B (orange). 




Figure 6: The dimcnsionless chiral condensate c versus the mass parameter m as defined 
in equation ([2.45)1 at baryon chemical potential \i B /M q = 0.01 (a), /j, B /M q = 0.1 (b) and 
[i B /M q = 0.8 (c) for the case Nf — 2. The dotted purple curve corresponds to Minkowski 
embeddings and the five other curves to black hole embeddings with isospin chemical 
potential [i 1 = (black), fi 1 = jfi B (green), — \pP (blue), fi 1 = jfi B (red) and 
/i 1 = \i B (orange). 




Figure 7: The dimensionless baryon density d versus the mass parameter m as defined 
in equation (|2.45| at baryon chemical potential fx B /M q — 0.01 (a), fx B /M q = 0.1 (b) and 
H B /M q = 0.8 (c) for the case Nf = 2. The five different curves correspond to black hole 
embeddings with isospin chemical potential fi 1 — (black), /j, 1 = jfi B (green), \i l = 
(blue), [i 1 = j[i B (red) and [i 1 = yP (orange). In the Minkowski phase, the baryon 
density is always zero and is therefore not shown in the figures. 
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(a) (b) (c) 

Figure 8: The dimensionless isospin density d 1 versus the mass parameter m as defined 
in equation (|2.45p at baryon chemical potential /i B /M q = 0.01 (a), /i B /M q = 0.1 (b) 
and ji B /M q = 0.8 (c) for the case Nf — 2. The five different curves correspond to 
black hole embeddings with isospin chemical potential [i 1 — (black), /i 7 = \\x B (green), 
[i 1 = (blue), /i 1 = j/j, b (red) and /i 7 = [l b (orange). In the Minkowski phase, the 
isospin density is always zero and is therefore not shown in the figures. For a zero isospin 
chemical potential /i 1 = 0, the isospin density d 1 is also zero and therefore coincides with 
the m-axis. 



embeddings at a finite mass m, the divergent behavior of the chiral condensate 
for large chemical potentials (/x B + p I )/Mq > 1 indicates that there is no phase 
transition between black hole and Minkowski embeddings anymore. 

Fig. [3 and [8] show the baryon and isospin density versus the mass parameter 
m calculated by the inversion of equation (|2,37[) at different baryon and isospin 
chemical potentials. These figures show that for all chemical potentials [i B l l , the 
baryon and the isospin density are zero at zero mass m = 0. For small chemical 
potential, fig. [3(a), (b) and El (a), (b) show that the densities increase until they 
reach a maximum at m ~ 1 . Beyond the maximum the densities decreases rapidly, 
but they do not reach to zero at the largest m which can be constructed by black 
hole embeddings. Therefore, there is also a discontinuous step in the densities as in 
the chiral condensate at the transition from black hole embeddings to Minkowski 
embeddings. For large chemical potentials (fi B + p, 1 ) / M q > 1 (fig. [3(c) andEl(c)), 
the densities diverge as we increase m, as does the chiral condensate. 

3 Canonical D7-brane thermodynamics 

3.1 The canonical ensemble: 

free energy, entropy, energy and speed of sound 

In this section we study the thermodynamic behavior of the quark sector at finite 
baryon and isospin density which is dual to thermal contributions of the D7-branes. 
Since the partition function Z of the field theory is given according the AdS/CFT 
dictionary by 

Z = e ~ Son - shcU , (3.1) 

with the Euclidean on-shell supergravity action Son-shell , the thermodynamical 
potential, in the canonical ensemble the free energy F, is proportional to the 
Euclidean on-shell action 

F = -T In Z = TS on . shcll . (3.2) 
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To calculate the thermal contributions of the D7-branes F7, we have to determine 
the Euclidean DBI action (|2.23p on-shell. 

First we perform a Wick rotation in the time direction to obtain the Euclidean 
DBI-actiorJl. Next we must renormalize this action by adding the appropriate 
counterterms I c t (see [52] for a review) since it diverges on-shell. Since these 
counterterms do not depend on the finite densities, we can write them as in [25,27], 

Jet = — ^ (pLx -m 2 ) 2 - Amc , (3.3) 

where /0 max is the UV-cutoff and 

T D7 V 3 vo\{&)N fe * H XN c N f V 3 T 3 
A/-a = ^ = 35 , (3-4) 

whit the Minkowski space volume V3. Then the renormalized Euclidean on-shell 
action Ir may simply be written as 

jj- = j^G(m, fi) - 1 ( / 9 min -m 2 ) 2 -4mc , (3.5) 

where p m i n determines the minimal value of the coordinate p on the D7-branes, 
i- e. p m \ n = 1 for black hole embeddings and p m m = Lq for Minkowski embeddings 
and 

G(m,p)= r dp lp 3 ff(l- X 2 ) 

J Pmin \ 



£ \ 1 " X 2 + P 2 (d P x) 2 - 2^(1 - X 2 )(^) 2 (3.6) 



x 

i=i 



f 

-Nfip^-pm 2 )^. 

Since this action depends on the background fields \ an d Xi and therefore 
on the mass parameter m and the chemical potential pi, it is proportional to the 
grand potential in the grand canonical ensemble. To get the free energy in the 
canonical ensemble, we must perform a Legendre transformation of the on-shell 
action Ir as we did in equation (|2.27p . The Legendre transformed on-shell action 
Ir is given by 

jj- = Jj-G(m, d) - i (p min - m 2 ) 2 - imc , (3.7) 

where 



G(m,d) 



r dp (p 3 ff (1 - x 2 ) ^l~X 2 + P 2 {d P xf 

J Pmin \ 



y 1 + N f U-pm 2 ) 



(3.8 



X 

i=l 



6 The Euclidean time must be periodic with period f3 = 1/T, such that the geometry is non- 
singular. 
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In the following we use this action to calculate the main thermodynamical quan- 
tities. The difference to the results obtained in [27] is induced by the sum over 
the Nf different densities in equation (|3.8|) . 



3.1.1 The free energy 

The contribution of the D7-branes to the Helmholtz free energy is 

F 7 = TI R = ^ F 7 (m, d) (3.9) 

with the dimensionless quantity 

F 7 (m,d) = = ^-G(m,d) - ^{(p 2 min - m 2 ) 2 - Amc} . (3.10) 

The dimensionless quantity J- 7 determines the dependence of the free energy F of 
the complete setup on the quark mass M q and the quark densities n q at a fixed 
temperature T. 

Using the example of the free energy we study how a thermodynamic quantity 
of the complete setup is composed of the quantities of the subsystems. The parts 
from the D3-branes F%, which is the known quantity for the conformal J\f = 4 
SYM theory, and from D7-branes Fj form the free energy of the complete setup 

F = F3 + F7 = 4^ + 3) 

Here we only get the first order contribution in Nf/N c since we work in the probe 
brane limit Nf <C N c . In the following we define quantities similar to T 7 for the 
other thermodynamic quantities and denote them also with a calligraphic letter. 
Note that these thermal contributions of the D7-branes vanish in the limit where 
supersymmetry and conformal invariance is restored. This limit is given by zero 
densities and large mass parameters m — > oo. 

3.1.2 The entropy 

As in [25,27], the contribution to the entropy by the D7-branes is given by 

dF 7 2 dF 7 XN c N f V 3 T 3 ~ 
Sr = ~ W = -*R ^ = 32 S 7 (m,d) , (3.12) 

with the mass and density dependent part 

Nt 

-12 , 

S 7 (m, d) = —4F 7 (m, d) + — - > difli — 2mc . (3.13) 

Nt — ' 

J i=l 

In comparison to the case without isospin density [27], there is sum over the Nf 
different densities. 
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3.1.3 The energy 

Using the thermodynamic relation E = F + TS, we calculate the energy 

E 7 = \ f 2 3 £r(m,d), (3.14) 
with the dimensionless quantity 

Nt 
12 ~ 

£7(171, d) = — 3^7(171, d) + — — y difli — 2mc . (3.15) 

J i=i 

Again a sum over the Nf densities appears in the equation above, which is due to 
the finite isospin densities. 

3.1.4 The speed of sound 

The speed of sound is given by 

2 _v 3 dP_ d(n 3 + n 7 ) or 

Vs ~ 8E ~ dT '^(E 3 + E7) , [ ' 

where the pressure P is given by PV 3 = — = —(^3 + Q7). Q and E are the 
grand potential and energy of the complete system. and E 3 are the contribu- 
tion to the grand potential and energy of the D3-branes, respectively. Using the 
thermodynamic relations S = — (^) VN = — (§P)y, an d c v = {^^vn 1 we ma ^ 
rewrite the speed of sound in the canonical ensemble 

2 S 3 + S 7 
v s = . (3.17) 

C-vZ > Cv7 

Using equation (|3.15|) the specific heat c V 7 becomes 

3XN c N f V 3 T 3 ( 1 2 dc 2^4-; u,r, n 

c, 7 = 35 7 \ mc - m — + — \d l [m—-^ l ) , (3.18) 



16 S -*» ■ W/ ^-V-* n 

J 1 = 1 

where we used 8t = —^d m (see equation (|2.45p ). With the contributions of the 
D3-branes given by the thermodynamics of the conformal M = 4 SYM theory, 

S3 = y V 3 N*T 3 and c v3 = 3S 3 , (3.19) 
we evaluate the speed of sound to first order in Nf/N c 



2 1 
vz = - 
s 3 



XN f 



(3.20) 



with 

Nt 

1 dc 2 ( dli- \ 

V S V, d) =rnc--m*- + £^ m J£-4*J . (3.21) 

Notice that the first order correction in Nf /N c is consistent with the approximation 
of probe branes (cf. with equation (|3.11|) ). The last term reflects the presence of 
the finite baryon and isospin density. The dependence on the mass parameter m 
and the chiral condensate c is given by the result without densities [25]. 
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(a) (b) 

Figure 9: The dimensionless free energy versus the mass parameter m as defined in 
equation (|2.45|) at baryon density d B = 5 ■ 10~ 5 (left) and d B = 2 (right) for the case 
Nf = 2. The five different curves in each figure correspond to d 1 = (black), d 1 = jd B 
(green), d 1 = \d B (blue), d 1 = § d B (red) and d 1 = d B (orange). 




m 



Figure 10: The dimensionless entropy Sj versus the mass parameter m as defined in 
equation (|2.45p at baryon density d B = 0.5 (left) and d B = 2 (right) for the case Nf = 2. 
The five different curves in each figure correspond to d 1 — (black), d 1 = jd B (green), 
d 1 = \d B (blue), d 1 = \d B (red) and d 1 = d B (orange). 

3.1.5 Numerical results for the free energy, entropy, energy and the 
speed of sound 

In fig. [TUl [TT1 and [12] we present numerical results for the mass and density 
dependent part of the thermodynamic quantities at different baryon and isospin 
densities for Nf = 2. 

Fig. [9] shows the free energy versus the mass parameter m at different den- 
sities as calculated from equation (|3.10j) . For small densities (see fig. 0(a)) we see 
again that the low temperature black hole embeddings have the same free energy 
as the zero-density Minkowski embeddings. Since the free energy determines the 
thermodynamics entirely, the same effect is also seen in the other thermodynamic 
quantities and we therefore do not consider them at small densities. For larger 
densities (see fig. El(b)) we see a displacement due to the isospin densities, which is 
constant over a wide range. Moreover, the free energy behaves linearly for m > 2. 
The increase of the free energy as m — > 00 demonstrates the deviation from the 
conformal M = 4 thermodynamics which was discussed below equation (|3.1ip . 
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(a) (b) 

Figure 11: The dimensionless energy £7 versus the mass parameter m as defined in 
equation (|2.45p at baryon density d B = 0.5 (left) and d B = 2 (right) for the case Nf = 2. 
The five different curves in each figure correspond to d 1 = (black), d 1 = jd B (green), 
d 1 = \d B (blue), d 1 = \d B (red) and d 1 = d B (orange). 




(a) (b) (c) 

Figure 12: The dimensionless speed of sound relative to the conformal case Vg versus 
the mass parameter m as defined in equation (|2.45[) at baryon density d B = 5 • 10 -5 (a), 
d B = 0.5 (b) and d B — 2 (c) for the case Nf = 2. The five different curves in each figure 
correspond to d 1 = (black), d 1 = \d B (green), d 1 = \d B (blue), d 1 = |J S (red) and 
d 1 = d B (orange). 

For small m, i. e. small quark masses compared to the temperature M q <C T, we 
see deviations from the linear behavior due to thermal fluctuations which increase 
the free energy. 

In fig. [10] the entropy £7 versus the mass parameter m as calculated from 
equation (|3.13p is shown. Again a displacement due to the isospin densities appear. 
However, the displacement in the entropy becomes smaller as we increase m and 
the entropy approaches a non-zero value independent of the isospin density as 
m — > 00. Moreover if the baryon density is above its critical value, a minimum 
in the entropy appears for zero isospin density. By increasing the isospin density 
this minimum disappears slowly. 

The energy £1 shown in fig.[TT]as calculated from equation (13.150 shows a very 
similar behavior as the free energy, discussed above. However, a minimum appears 
in the energy at m ~ 1.3 for a baryon density above the critical density (see fig. [TTl 
(a)). Increasing the baryon densities (see fig. [TTJ (b)), this minimum disappears. 

Fig. [I2l shows the speed of sound Vf calculated from equation (I3.2ip . Again 
a displacement due to the isospin densities appear. In the speed of sound this 
displacement grow as we increase m. For a baryon density above the critical 
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density and zero isospin density, the speed of sound is constant at m « 1.3 as 
shown in fig. [12] (b). Increasing the baryon or isospin density the speed of sound 
monotonously decreases as m increases (fig. E] (b) and (c)). 

3.1.6 Interpretation: Quark mass selection and temperature change 

Here we discuss the effect of a non-zero isospin density on the quark mass M q . We 
argue that under certain conditions explained below, the thermodynamic quanti- 
ties given in the previous section 13.1.51 select preferred values M* for the mass M q 
of the quarks inside the quark gluon plasma. In a slightly different setup we find 
that our system prefers a distinct temperature T = T*. 

In our model the quark mass M q is a free parameter. We fix the tempera- 
ture T and charge densities n l and keep the volume V of our system constant. The 
canonical ensemble describes this situation with the free energy being the relevant 
thermodynamic potential. We argue how the quark mass M q as a free parameter 
has to be selected in order to satisfy the second law of thermodynamics. In fact 
we will find a competing behavior between entropy and energy requirements. The 
second law of thermodynamics implies that in thermal equilibrium, the entropy is 
maximized at fixed thermodynamic potential. This is equivalent to minimzing the 
thermodynamic potential at a fixed entropy [53]. In our case this potential is the 
free energy. However, competing situations between entropy and energy require- 
ments are well known from frustrated systems, e. g. polymers or anti-ferromagnets 
on a triangle lattice, as well as glass states, e.g. spin glass [54]. In these systems 
the temperature is usually the thermodynamic variable which adjusts at a distinct 
value where the criteria are optimally satisfied. In our setup we have the quark 
mass M q as a free parameter, which in analogy to the case discussed above adjusts 
at a specific value M*. In order to find this optimal value M*, we examine fig- 
ures M and [TUl Recall from equations (|2.45p and (|2.35j) that the mass parameter m 
is related to the temperature and quark mass by m = 2M 9 /(VAT) and n l oc d l . 
We observe that in the region m > 4 for all densities, the free energy decreases at 
a constant entropy as the quark mass M q is increased. At smaller values of the 
mass parameter m < 4, the entropy changes as we vary the quark mass M q . The 
second law of thermodynamics forbids the system to be in equilibrium. Prom these 
considerations and from looking at figures [9] and [10] we conclude that the quark 
mass has to take a value of < M* < \/ATto/2, such that < m < 4. Since in 
this regime the entropy and free energy both vary when the densities are changed, 
we expect the value of M q to depend on the densities n\ i.e. M q = M*(n l ). There- 
fore, our system with the quark mass as a free parameter shows similar properties 
as the frustrated systems discussed above. 

It is amusing to apply the results of this analysis to the quark gluon plasma 
formed after the big bang. This would imply that our universe could have evolved 
in such a way that the quark mass was chosen such that it optimizes the free 
energy and entropy requirements. Note that during this mass selection process we 
have to stay in a regime where the metric does not change considerably in order 
for this system to be approximately closed [53]. 

Now we consider a different scenario, by which a distinct temperature T* is 
favored. We work in the microcanonic ensemble by fixing the entropy instead of 
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Figure 13: (a) The dimcnsionless free energy ^7, (b) the speed of sound V, , (c) the chiral 
condensate c, (d) the baryon chemical potential yP /M q and (e) the isospin chemical 
potential [i 1 /M q versus the mass parameter m as defined in equation (|2.45p near the 
phase transition at the baryon density d B = 5 • 1CP 5 for Nj = 2. The resolution of 
the figures is not high enough to dissolve the curves at different isospin density d 1 = 
0, \d B , \d B , jd B , d B . The dashed line at m = 1.306 marks the phase transition. 



temperature. An example of such a situation is the collision of heavy ions in vac- 
uum. Since the natural temperature variable for our system is the temperature T 
measured in units of the quark mass M q , we let m change freely. In this case 
the quark mass M q is fixed and the volume of plasma evolves to a preferred value 
of temperature T, in a reversible thermodynamic process. The system seeks to 
minimize its internal energy £7 while its entropy £7 is fixed by definition. In this 
setup the plasma relaxes isentropically to a favored temperature T = T* . 

3.2 Phase transition/phase diagram in the canonical ensemble 

In the pure baryonic case [27] it was shown that there is a first order phase tran- 
sition between two black hole embeddings for densities below the critical density 
^crit = 0.00630 In this section we analyze how the isospin density influences this 
phase transition. In fig. [13] we give an example of the behavior of the free energy, 
the speed of sound, the chiral condensate, the baryon and the isospin chemical po- 
tential [i B l l /M q close to the phase transition. The free energy, the speed of sound 
and the chiral condensate behave as in the pure baryonic case [27]. The isospin 
chemical potential fi 1 /M q takes different values depending on the isospin density 
after the phase transition (see fig. [13] (e)). The baryon chemical potential (fig. [131 
(d)) instead is almost independent of the isospin density. However, the behavior of 
both chemical potentials changes dramatically if the baryon and isospin densities 
are equal in agreement with equation (12.420 . 

Fig. [13] shows that the phase transition is first order since the free energy 

7 Note that the normalization we use in this paper (|2.28jl differ to the normalization used 
in [27] by N f . 
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Figure 14: (a) The phase diagram in the canonical ensemble for Nf = 2: The cusp 
marks the phase transition for d B — d 1 — and the upper rim the critical points. Notice 
that the phase diagram is approximately 0(2) invariant in agreement with the symme- 
tries discussed in section 12.3.11 (b) Position of the critical points: Comparison of the 
position of the critical points (cross) with the 0(2) invariant circle (blue) and a superel- 
lipse (<F) a = (df rit ) a - {d B ) a with a = 2.08 (dashed). The difference of the parameter a 
from 2 determines the deviation from the 0(2) symmetry. 

is continuous and the speed of sound, the chiral condensate and the chemical 
potentials show a multivalued behavior and therefore a discontinuous step at the 
phase transition. The phase transition is marked by the crossing point in the free 
energy of the two branches coming in from small m and large m (see fig. [T3] (a)). 

In fig. (a) we map out the phase diagram by varying the densities and 
determining the crossing point in the free energy. As outlined in section 12.3.11 the 
phase diagram is approximately 0(2) invariant. We find a line of critical points, 
which mark the critical mass parameter m cvl t and critical densities d^ it , d^ lit . 
For densities around the critical points, we expect the largest deviation to the 
0(2) symmetry in the phase diagram (see equation (|2.30p ). To get a quantitative 
description of this deviation, we plot the critical points and compare them to the 
0(2) invariant circle and the superellipse (d 1 ) 0, = {d B it ) a — (d B ) a with a = 2.08 in 
fig- ttH(b). The deviation a — 2 = 0.08 is still small at the critical points. Thus, the 
complete phase boundary in the canonical ensemble is 0(2) invariant in a good 
approximation . 

More advanced thermodynamical investigations of the pure baryonic case [27, 
28] show that there is a region were the black hole embeddings are not a stable 
ground state for this theory. The points where the phase transition between the 
two black hole embeddings occurs is at the boundary of this unstable region. 
In [28] the authors expect that in this region an inhomogeneous mixture of black 
hole and Minkowski embeddings is the stable ground state. Therefore the phase 
transition is not between two black hole embeddings, but between a black hole 
embedding and this inhomogeneous mixture. The existence of a mixed phase is 
typical for a first order phase transition and can e. g. be observed in boiling water. 
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(a) (b) 




(c) (d) 

Figure 15: (a) The baryon chemical potential fi B /M q , (b) the isospin chemical potential 
pj /M q , (c) the chiral condensate c and (d) the free energy T7 versus the baryon and 
isospin density at m = 20, where the mass parameter m is defined in equation (|2.45[) . 
in the case Nf = 2. (c) and (d) are contour plots of the chiral condensate and the free 
energy, respectively. 

The water does not instantly turn into gas but forms droplets consisting of a 
mixture of water and water vapor. 

3.3 Canonical thermodynamics in the large quark mass limit 

In the large mass limit m — * 00, i. e. M q 3> T, the dimensionless thermodynamic 
quantities and the chiral condensate vanish if the baryon and isospin densities are 
zero because supersymmetry and conformal symmetry are restored and therefore 
the thermodynamics are given by the results of the stack of D3-branes dual to the 
conformal N = 4 SYM theory [25]. Introducing finite baryon and isospin densi- 
ties break these symmetries even though the AdS black hole background becomes 
symmetric and we get deviations in the thermodynamic quantities from the con- 
formal results (see e. g. fig. [9]). In the following we investigate the dependence of 
these deviations on the densities d * further. 
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Figure 16: The phase diagram for m — > oo, i.e. M q ^> T: The dashed lines marks the 
phase transitions where the chemical potentials develop a discontinuous gap. Notice that 
the dashed lines are the fixed points of the symmetries d\ <-*■ —d\, di <-> — c?2- 



3.3.1 Numerical results 

First we discuss the numerical results for the thermodynamic quantities, the chem- 
ical potentials and the chiral condensate at a fixed mass m = 20 in dependence 
of the baryon and isospin density for Nf = 2, presented in fig. [151 We restrict 
our analysis to non-negative baryon and isospin densities dP^d 1 > 0. The results 
may be continued into the other three quadrants by the accidental symmetries, 
discussed in section [2.3.1i The discontinuous step in the chemical potentials along 
the line d 1 = d B > described in and (JX25|) are shown in fig. [TS] (a) and 

(b). We can easily characterize this behavior in the large mass limit m — > oo by 

!L. ^H(P-dF) and iL _> ff(( p _ JB) , (3.22) 

where H is the Heaviside function. This discontinuity suggests that there is a first 
order phase transition between the two regions < d 1 < d B and < d B < d 1 . At 
the same line the chiral condensate and the thermodynamic quantities develop a 
kink (see fig. [T5j) , which also indicates that there is a phase boundary along the line 
d B = d > 0. Using the accidental symmetries discussed in section [2.3.11 we give 
an overview of the four different phases and the values of the chemical potentials 
in the corresponding region in fig. [T6l The same four different phases are also seen 
in QCD-like theories studied in [46]. Since in [46] the authors study the grand 
canonical ensemble, the phase transition is marked by a discontinuous step in the 
baryon and isospin density and phase boundaries are given by \n B \ = Along 
these phase boundaries the chiral condensate develops a similar kink as we observe 
here. 

From the contour plot in fig. [T^] (c), we see that the chiral condensate is 
independent of the isospin density for large masses in the region < d 1 < d B 
since the contour lines are parallel to the lines of constant baryon density in 
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the region < d 1 < d B . To investigate the region < d B < d 1 , we can simply 
interchange the baryon and isospin density by the accidental symmetries discussed 
in section 12.3.11 Therefore, the chiral condensate only depends on the isospin 
density in the region < d B < d 1 . Fig. [3] shows that the chiral condensate is 
also independent of the mass parameter m since it approaches a constant value as 
m — > oo. Approximately, the dependence of the chiral condensate is given by 

for < d 1 < d B 

as to — > oo . (3.23) 

for < d B < d 1 

Next we analyze the thermodynamic quantities and start with the free energy. 
The contour plot of the free energy in fig. [JS] (d) directly shows a dependence on 
the isospin and the baryon density in both regions < d 1 < d B and < d B < d 1 . 
In the region < d 1 < d B , the dependence appears to be linear on the baryon 
density but nonlinear on the isospin density. Fig. [9] shows also a linear dependence 
of the free energy on the mass parameter to in the large mass limit. In order collect 
this behavior in a formula, we make a suitable ansatz and fit to the numerical data 

( V2md B - 2d B + 0.6§^J for < d 1 < d B 
^7~< „ „ ji 2 2 5 ~ ~ as m^oo. (3.24) 

[ V^md 1 - 2d 1 + 0.6 y 9 for < d B < d 1 

Similarly we write for the energy 

( V2md B + 0.6^p^5 for < d 1 < d B 
£7 ~ \ ~ i,ib\2 25 ~ ~ as to — > 00 . (3.25) 

{ V^md 1 + 0.6iij^ for < d B < d 1 

The entropy behaves similarly as the chiral condensate and we therefore may 
approximate it with 

f 2d B for < d l < d B , 
07 — > < ~ T ~„ ~ T as to — ► 00 . (3.26) 

[2d 1 for < d B < d 1 

In the following we interpret the large mass limit of the chiral condensate, the free 
energy, the energy and the entropy found in (|3.24|) - (|3.26|) . 



3.3.2 Interpretation of the numerical results 

For the interpretation of the results found in (|3.24|) - (j3.26|) . the difference of the 
two regions < d 1 < d B and < d B < d 1 as displayed in fig. [15] is essential. First 
we show that in both regions the larger density is a measure for the degrees of 
freedom. For this purpose, we map our theory to a QCD-like theory. In our setup 
the open D3-D7 strings model quark^l of a QCD-like theory [22]. For Nf = 2 we 
consider two different quarks, e. g. u and d quarks. We assign a baryon charge q B 
and isospin charge q 1 to these quarks, 

u = \q B ,q I ) = \1,1), d=|l,-l). (3.27) 
8 Strictly speaking the D3-D7 strings model a SUSY multiplet consisting of quarks and squarks. 



31 



The effect of the symmetries (I2.33D on these charges is given by 
d% <-> d,2 ■ 

u= |1,1) i-> |1,-1) = d and d = |1, -1} ■-»• |1, 1) = u 

d\ <-» — di : 

u = |1, 1) i— > | — 1, —1) = u and d = |1, —1) i— > |1, —1) = d 
d 2 ^ - d 2 : 

u = |1, 1) i— > |1, 1) = u and d = |1, — 1) i— >• | — 1, 1) = <i , 

(3.28) 

where we defined the anti-quarks u and d consistently with the ideas given in 
section [2.3.11 and the convention used in QCD-like theorie^. Using the charges of 
the quarks we may assign a baryon charge to the plasma by 

plasma = n u + n d - n u - n d , (3.29) 

and the isospin charge by 

<?plasma = n u + % ~ Tift - Tl d , (3.30) 

where rii counts the different quarks in the plasma. A plasma with < f/pi asma < 
plasma' < d 1 < d B , may be constructed using u and d quarks only. Ex- 
perimentally such a plasma can be produced by heavy ion collisions. However, 
for a plasma with < qp[ asma < 9pi asma ) i- e. < d B < d 1 , we must consider the 
free quarks u and anti-quarks d. Therefore for positive densities, the larger one of 
the two quantities </p^ asma and <7p lasma counts the degrees of freedom of the plasma 
consisting of the free quarks and anti-quarks. This explains the difference in the 
two regions displayed in fig. [T5l 

First we consider the entropy 1S7 given in equation (|3.26p . The entropy mea- 
sures the logarithm of the number of states the plasma could be in. The degrees 
of freedom increase these number of states by their phase space volume. Since 
the larger density determines the degrees of freedom, it is clear that the entropy 
increases linearly with the larger density. 

Next we consider the free energy and the energy found in (|3.24|) and (|3.25[) . 
Since the degrees of freedom contribute to the free energy (energy) by their 
rest mass, a term appears which is proportional to d B m and d l m. In the free 
energy (energy) there is also a nonlinear dependence on the smaller density, 
(d 1 ) 2 - 25 /(d B ) - 9 and (d B ) 2 - 25 /(d 1 ) - 9 , respectively. We expect that this contribu- 
tion is due to a charge of the plasma which is induced by the non-Abelian charges 
of the quarks and measured by the smaller density. This charge is similar to an 
electrostatic charge in the sense that the free energy (energy) increases as we in- 
crease the amount of the charge of the plasma. The deviation in the exponent of 
(d 1 ^) 2,25 from the naively expected Coulomb behavior (d J / B ) 2 is probably due 
to non-perturbative effects of vacuum polarisations since the gauge coupling is 
largj^l. The fact that there is a term depending on the larger density (J B )~ ' 9 
may be due to a screening of the plasma charge by the quarks and anti-quarks. 

9 Recall that the pion ir + — ud — |0, 2} (tv~ = du — |0, —2)) has positive (negative) isospin 
charge in QCD. Notice that the normalization of our charges is not the same as in QCD. 

10 In QED, where the coupling constant is small, the perturbative corrections due to vacuum 
polarisations are sometimes called Uehling effect [55]. 
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3.3.3 Expected behavior of the system perturbed by interactions 
which break baryon and isospin conservation 

In this section we study the possible behavior of our system if the baryon and 
isospin density are not conserved. Interactions which break baryon (isospin) con- 
servation are e.g. included in the MSSM (electroweak theory). We are not able 
to include such interactions explicitly in our setup. However, we may analyze the 
modification in the plasma induced by variation of the density using the thermo- 
dynamic results found in (|3.24p to (|3.26p . The variation in the plasma must be 
such that the free energy becomes minimal at a fixed entropy (see section [3.1. 6ft . 
We discuss the reaction of our system separately for the four different regions in 
the phase diagram shown in fig. [TBI 

Experimentally the region d B > \d?\ may be observed by heavy ion collisions. 
In this region we learn from equation (|3.26p that the entropy is independent of 
the isospin density but proportional to the baryon density. Therefore, at constant 
entropy the baryon density must be constant, too. Only the isospin density may 
change. For constant baryon density the minimum of the free energy is located 
at zero isospin density, see fig. [15] (d). Thus, we expect our system, perturbed by 
weak interactions, to decrease the amount of isospin density to zero. This effect is 
well known in nuclear physics. The nuclei stable with respect to /3-decay consists 
of the same number of protons and neutron^]. As is well known the /3-decay is 
induced by electroweak interactions, which we consider here to change the isospin 
density. 

The region d B < — Id 1 ] is similar to the region discussed above since they 
are connected by the symmetry d B <-» — d B which interchanges particles with 
anti-particles. Thus, this region determines the physics in a situation where anti- 
particles dominate. In this region we also expect that the system approaches 
a state of zero isospin density. For the example of nuclear physics this means 
that nuclei are also stable under (3 decay if they consist of the same number of 
anti-protons and anti- neutrons. 

In the other two regions the amount of isospin charges is bigger than the 
amount of baryon charges \d B \ < Id 1 ]. To relate these two regions with the regions 
discussed above we make use of the symmetry d B <-» d 1 which interchanges the 
baryon and isospin densities. Thus, we expect that for these regions the system 
moves to a state of zero baryon density at a fixed isospin density. Such a state can 
be achieved by pairs of particle and anti-particle of the same isospin since these 
pairs only contribute to the isospin density but not to the baryon density. An 
example for such pairs are mesons as e. g. ir + and tt~ . 

3.4 Canonical thermodynamics in the small quark mass limit 

In the small mass limit m — > 0, i. e. M q <C T, the thermal energy is sufficient 
to excite the plasma, e. g. by producing pairs of particles and anti-particles. The 
probability of these thermal excitations is in the large temperature limit, i. e. small 
mass limit, approximately described by the Boltzmann factor 
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For heavy nuclei there is a deviation due to the electric charge of the protons. 



33 




Figure 17: The deviation in the free energy due to thermal fluctuations ^t hormal versus 
the mass parameter m as defined in equation (|2.45p at a baryon d B = 0.5 (left), d B = 2 
(right) and different isospin density d 1 for Nf = 2 in a logarithmic plot. The different 
colors correspond to different isospin densities d 1 — (black), d 1 = \d B (green), d 1 = \d B 
(blue), d 1 = jd B (red) and d 1 = d B (orange). 



where E is the energy of the plasma. In our setup this expression may be approx- 
imated by 

P = e -Am+BY,id I if l h( m ) ^ (3.32) 

where A, B are some constants. The dependence of the chemical potentials fx 1 * on 
the mass parameter m induces a slow variation of 

K{m) = B^2 d h £l u (m) , (3.33) 

i 

since the dimensionless chemical potentials fl 1 " 1 are almost constant for small 
masses m [27, fig. 8(a)]. Thus, the Boltzmann factor may be written as 

p = e -Am+K(m) ? (3 34) 

where K(m) is slowly varying. We expect that the deviation in the thermodynamic 
quantities due to thermal fluctuations, i. e. the deviation from the large mass limit 
studied in section 13.3} behaves as the Boltzmann factor described in equation 
(|3,34p . In the example of the free energy the thermal fluctuations ^ hermal are 
given by the dimensionless free energy calculated in equation (|3.10p minus the 
large mass limit of the free energy given in equation (I3.24|) . 

f V2md B - 2d B + 0.6^JJ for d B > d 1 

^-thermal T ) (" ) ' (o or:\ 

( V^md 1 - 2d 1 + 0.6 ( y )0 .9 for d < d 

A selection of our numerical results at different baryon and isospin densities is 
given in the logarithmic plot in fig. [T71 where we see the expected result at small 
masses m. We also observe deviations in the large mass limit, i. e. T <C M q , since 
the Boltzmann factor given in equation (|3.3ip does not give the full quantum 
statistic needed to describe the system. 
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4 Grand canonical D7-brane thermodynamics 



4.1 The grand canonical ensemble: 

Grand potential, entropy, energy and speed of sound 

In this chapter we investigate the properties of the grand canonical ensemble. 
For this purpose, we have to use the renormalized Euclidean action (|3.5[) . It is 
proportional to the grand potential 

\N r NfViT 4 
V 7 = TI R = 32 W ^ A) . 

with the dimensionless quantity 

Wr(m, ft) = j^G(m, Jm) - - [(^ in - m 2 f - 4mc] . (4.2) 

Since we must also consider Minkowski embeddings in the grand canonical ensem- 
ble, we rewrite G(m, fj) from equation (13. 6p in the coordinates L and r suitable 
for Minkowski embeddings and for constant gauge fields 

/•oo r , 

G{m) =N f dr ffr z ^l + (d r L) 2 + (r + Ld r L) [m 2 - (r 2 + L 2 )] . (4.3) 
Jo 1 J 

For Minkowski embeddings, the grand potential does not depend on the chemical 
potential and therefore coincides with the free energy at zero chemical potential 
[25]. Since the free energy determines the thermodynamics entirely, we may apply 
the results of [25]. 

We calculate the thermodynamic quantities in the black hole phase numerically 
by evaluating the canonical results at a fixed chemical potentials. The transfor- 
mation to fixed chemical potential can again be done, as in section 12.51 by the 
functions di(p,i,m), which is the inverse of (12.371) . Only the speed of sound has 
to be recalculated as it is no thermodynamic state function. The difference to 
the result (I3.2ip in the canonical ensemble is that the derivative now acts on the 
density instead of the chemical potentials, 
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V s (m, M ) = mc - -m — + — ^ H (m— - Ad, j (4.4) 

since the specific heat in the grand canonical ensemble is determined by cy = 
iw)v)j, ( c ^' derivation in section I3.1.4|) . 

4.1.1 Numerical results for the grand potential, entropy, energy and 
speed of sound 

In fig. [TSl [T9l and [20] we present numerical results for the thermodynamic quan- 
tities in the grand canonical ensemble for Nf = 2 at different baryon and isospin 
chemical potentials. 

Fig. [18] shows the dimensionless grand potential W7 versus the mass parameter 
m as calculated by equation (|4.2|) . we compare the grand potential given in fig. [TH1 
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(a) (b) (c) 



Figure 18: The dimensionless grand potential WV versus the mass parameter m as 
defined in equation (|2.45j) at baryon chemical potential p B jM q = 0.01 (a), pP jM q = 0.1 
(b) and p B /M q = 0.8 (c) for the case Nj = 2. The dotted purple curve corresponds to 
Minkowski embeddings and the five other curves to black hole embeddings with isospin 
chemical potential fi 1 = (black), // = \\x B (green), // = \p B (blue), p 1 = \p B (red) 
and [i 1 — ji B (orange). 
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(a) (b) 

Figure 19: The dimensionless entropy £7 versus the mass parameter m as defined in 
equation (|2.45p at baryon chemical potential pP jM q = 0.1 (left) and \i B jM q = 0.8 (right) 
for the case Nt = 2. The dotted purple curve corresponds to Minkowski embeddings and 
the five other curves to black hole embeddings with isospin chemical potential p 1 = 
(black), p 1 = jp B (green), p 1 = ^/i B (blue), p 1 = jp B (red) and p 1 = p B (orange). 

with the result where no finite charges are considered [25]. The deviation between 
these two cases is induced by the finite densities given in fig. \7\ and El Since the 
densities d 1 ' are all zero for m = 0, there is no deviation. By increasing the mass m, 
the densities and the deviation grow. At m ~ 1 the densities are maximal for small 
chemical potentials and therefore, the deviation becomes smaller as we increase 
m further and we get a phase transition similar to the case without chemical 
potentials (see fig.lTSI (a)). For larger chemical potentials but still (p B + [i I )/M q < 
1 (fig. [18] (b) and black curve in fig. [18] (c)), the deviation to the case without 
chemical potential is much larger and the position of the phase transition moves 
to larger mass m. For {ji B + ir)/M q = 1 (green curve in fig. [18] (c)), the grand 
potential approaches a constant value in the limit m — * 00. This value is smaller 
than the value in the Minkowski phase, such that there is no phase transition 
between black hole and Minkowski embeddings possible. For even larger chemical 
potentials (p B + [i I )/M q > 1, the grand potential diverges to —00 as m — > 00. We 
expect the black hole embeddings to be thermodynamically favored for all masses 
m if (fi B + p I )/M q > 1. 
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(a) (b) (c) 

Figure 20: The dimcnsionless speed of sound relative to the conformal case versus the 
mass parameter m as defined in equation (I2.45[) at baryon chemical potential (x B /M q = 
O.Of (a), [i B /M q = O.f (b) and [L B jM q = 0.8 (c) for the case N f = 2. The dotted 
purple curve corresponds to Minkowski embeddings and the five other curves to black 
hole embeddings with isospin chemical potential /j, 1 = (black), fj, 1 — jfi B (green), 
[i 1 = (blue), fi 1 = (red) and /j, 1 — [i B (orange). 

Fig. [19] shows the entropy S7 versus the mass parameter m. For small chemical 
potentials, the entropy behaves as the zero chemical potential result from [25]. By 
increasing the chemical potentials the spiral behavior disappears and the entropy 
decreases monotonically (see fig. [19] (a)). For large chemical potentials presented 
in fig. [19] (b) the entropy changes its behavior dramatically. It increases to a 
maximum and decreases again towards the values obtained for the Minkowski 
embeddings as we increase the mass m. For even larger chemical potentials the 
entropy diverges to 00 as m — > 00. The energy behaves like the entropy and we 
therefore do not present it separately. 

Fig. [20] shows the speed of sound Vg as calculated from equation (14. 4p . In the 
speed of sound there is no dramatic change aside from the transition at (fi B + 
[i^/Mq = 1 as we increase the chemical potential. By increasing the mass m it 
always decreases until the phase transition and increases again in the Minkowski 
phase. The magnitude of the speed of sound close to the phase transition increases 
as we increase the chemical potentials. 

4.2 Phase transition/phase diagram in the grand canonical en- 
semble 

We review the regions in the phase diagram for the pure baryonic case, where 
Minkowski and black hole embeddings can be constructed. From [23, 25] it is 
known that Minkowski embeddings can only exist for m > 1.3. In [28] it is shown 
that there is a region (large m and small fi/Mq), where black hole embeddings can- 
not exist (see also our numerical results in fig.[H(a) and (b)). Nevertheless, there 
is a region where Minkowski and black hole embeddings both can be constructed, 
but only one or a mixture of both is thermodynamically favored. In fig. [22] (a) we 
sketch these regions. 

In the following we determine the position of the phase transition numerically. 
For this purpose, we look for a crossing point of the grand canonical potential 
calculated using black hole and Minkowski embeddings. Fig. [2H shows the free 
energy, the speed of sound and the chiral condensate close to this transition. The 
entropy and the energy show a similar behavior as the chiral condensate. 
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(a) (b) (c) 

Figure 21: (a) The dimensionless grand potential Vty, (b) speed of sound V|f and (d) 
the chiral condensate c versus the mass parameter m as defined in equation (|2.45|) near 
the phase transition at fi B /M q = 0.01 for Nf = 2. The different colors of the curves 
corresponds to the different isospin chemical potentials [i 1 = (black), [i 1 = j/i b (green), 
[i 1 = \[i B (blue), /j, 1 = ^/jI B (red) and fi 1 = /i B (orange). The purple curve corresponds 
to Minkowski embeddings. 
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Figure 22: (a) Sketch of the possible embeddings in the grand canonical ensemble for one 
brane: In (A) only Minkowski embeddings, in (B) Minkowski and black hole embeddings 
and in (C) only black hole embeddings can exist. In (B) only one or a mixture of the two 
embeddings can be thermodynamically preferred. 

(b) Contour plot of the phase diagram in the grand canonical ensemble for Nf — 2: Con- 
tours label the value m _1 of the phase boundary. For small chemical potential the phase 
diagram is approximately 0(2) invariant in agreement with the symmetry discussed in sec- 
tion For larger chemical potentials the rotation symmetry breaks down and there 
is only an Z4 symmetry left. For (ii b + /j, I )/M q > 1, there is no transition between black 
hole and Minkowski embeddings. Below the hypersurface the Minkowski embeddings are 
thermodynamically favored. Above the hypersurface black hole embeddings are preferred. 
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By varying the chemical potentials, we map out the phase diagram for the 
grand canonical ensemble shown in Fig. [22] (b). For lower values of m" 1 , the 
Minkowski embeddings are the thermodynamically favored embeddings and for 
larger ones above the hypersurface the black hole embeddings are preferred. The 
contours at small chemical potentials show that the phase diagram is approxi- 
mately 0(2) invariant, in agreement with the discussion in section 12.3.11 This 
rotation symmetry breaks down for larger chemical potentials and the discrete 
Z4 symmetry is left. For (fi B + /i )/M q > 1 the black hole embeddings are the 
thermodynamically preferred embeddings for all masses m. Therefore, there is 
no phase transition between black hole and Minkowski embeddings for chemical 
potentials (fi B + fi I )/M q > 1. This is the natural extension of the pure baryonic 
result [28, 51] considering the sum of the baryon and isospin chemical potential 
fi B + (J, 1 as the relvant quantity. 

Our numerical results suggest that the phase transition between black hole 
and Minkowski embeddings is first order for finite chemical potentials as in the 
case of zero chemical potentials [23,25]. However, in [51] the authors show analyt- 
ically that the phase transition in the pure baryonic case is second order at zero 
temperature. We expect that the same is true if we include the isospin charges. 
Therefore, the order of the phase transition must change from first to second order 
at one point of the phase boundary. Due to numerical errors we cannot strictly 
distinguish between first and second order phase transitional. However, we argue 
in the following that this point should be located at zero temperature. In general 
at a critical point, where a phase boundary ends, a first order phase transition 
becomes second order [47]. We do not observe that the phase transition disap- 
pears. Thus, a critical point can only be located a the end of the phase boundary, 
i. e. at zero temperature or zero chemical potentials. Since the phase transition at 
zero chemical potentials is clearly first order [23,25], the critical point can only be 
located at zero temperature. Therefore, we expect a second order phase transition 
at zero temperature. To strengthen this argument one has to expand the analytic 
studies in [51] to the next order in the temperature T and determine analytically 
the order of the phase transition for a small but finite temperature T. 



4.3 Grand canonical thermodynamics in the large quark mass 
limit 

In the large quark mass limit m — ► 00, i.e. M q ^> T, the Minkowski embeddings 
are thermodynamically favored for small chemical potentials (fi B + n I )/M q < 1. 
Since in the Minkowski embedding the chiral condensate and the thermodynamic 
quantities are independent on the chemical potentials, they all vanish for m —* 00 
as in the case of zero chemical potentials [25]. The asymptotic behaviour is given 
by 

12,7 1 

yy 7 _ ; 57 ~ — j , £7 ~ j , c ~ =■ . (4-5) 

12m 4 3m 4 12m 4 6m & 

However for chemical potentials greater than the critical values, black hole em- 
beddings with finite baryon and isospin densities are preferred. From the studies 



It is numerically very challenging to construct all black hole embeddings close to the phase 
transition. 
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(a) (b) 

Figure 23: (a) The baryon density d B , (b) the isospin density d 1 versus the baryon and 
isospin chemical potential \i B /M g , ji 1 /M q at a fixed mass parameter m = 20, the mass 
parameter m is defined in equation (|2.45p . 




(a) (b) 

Figure 24: (a) The chiral condensate c and (d) the energy £7 versus the baryon and 
isospin chemical potential \i B /M q , /j, 1 /M q at a fixed mass parameter m = 20, where the 
mass parameter m is defined in equation (12.451) . 
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of the canonical ensemble, we know that the finite densities break supersymme- 
try and conformal invariance and therefore the thermodynamic quantities and the 
chiral condensate do not vanish in the large mass limit any more. 

In this section we study the system in the grand canonical ensemble and com- 
pare it with the canonical ensemble, which we studied in section [231 In fig- 1231 and 
Owe present the baryon and isospin density, the chiral condensate and energy at 
a fixed mass m = 20 versus the baryon and isospin chemical potential. The grand 
potential and the entropy behave like the energy or the chiral condensate. The 
phase transition between Minkowski and black hole embeddings is along the line 
{jjL B +[i I )/M q 1. Since the quantities do not depend on the chemical potentials in 
the Minkowski embeddings, we see a plateau in the region (fj, B +/j, I )/M q < 1 where 
all the presented quantities in fig. [23] and [23] are approximately zero. After the 
phase transition the baryon and isospin density increase rapidly (d B , d 1 ss 8 • 10 4 
at p B /M q = fi 1 /M g = 2.5 see fig. I23|) . Since we learn from the investigation of the 
canonical ensemble that the thermodynamic quantities and the chiral condensate 
depend at least linearly on the densities, the quantities also increase or decrease 
very fast after the phase transition (see fig. 124"]) . 

In the contour plots of the baryon and isospin density (see fig. [23]) we see 
a turning point in the contour lines at = y} for small chemical potentials. 
Since there is a phase transition between the regions < d 1 < d B and < 
d B < d in the canonical ensemble, we also expect a phase transition between the 
regions < yJ < y B and < fi B < fj, 1 in the grand canonical ensemble. In the 
canonical ensemble this phase transition is clearly marked by a discontinuous step 
in the chemical potentials, and by a kink in the thermodynamic quantities and 
the chiral condensate (see fig. [To]) . In the grand canonical ensemble this phase 
transition is not visible in the contour plots of the baryon and isospin density. 
The thermodynamic quantities and the chiral condensate also appear to behave 
smoothly at this point (see fig. !24p . Close to the line fi B = , the thermodynamic 
quantities depend only on the sum of baryon and isospin chemical potential, pP + 
p} . It appears that there is an inconsistency in the transformation between the 
grand canonical and canonical ensemble since we only observe a phase transition 
in the canonical ensemble. 

In the following we resolve the inconsistency and show that a new phase must 
be included in the grand canonical phase diagram (see fig. [22]) . For this purpose, 
we determine the different regions in the grand canonical phase diagram fig. [22] 
which the canonical and the grand canonical ensemble control. In the canonical 
ensemble, the possible values for the chemical potentials yPI 1 are 0, M q /2 and 
M q (see fig. dn] and equations (|2.41|) . (|2.42|) ) in the limit m — ► oo. Different 
states in the canonical ensemble determined by (d B ,d I ,m — * oo) are mapped to 
a single state in the grand canonical ensemble, e.g. (y B = Mq,^ 1 = 0, m — » oo). 
Therefore, the phase transition observed in the canonical ensemble (see fig. [To]) 
in the limit m — > oo cannot be seen in the grand canonical ensemble since the 
two dimensional density plane is mapped to three different points in the grand 
canonical ensemble (fi B = Mq,^ 1 = 0), (fi B = M q /2,p T = M q /2) and (p B = 
0,/j 1 = Mq). These three points are located along the phase boundary of the 
Minkowski phase (fi B + \i )/Mq ~ 1 in fig. [23] and [Ml However, we can analyze 
the plots in the grand canonical ensemble in this section for finite mass. We 
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also observe in the canonical ensemble that for finite mass and large densities 
the phase transition in the density plane disappears slowly, such that there is no 
inconsistency between the canonical and grand canonical ensemble in the region 
of finite mass m. 

Since we cannot map bijectively the canonical to grand canonical ensemble in 
the region of large mass m and chemical potentials pPI 1 ~ M q , the considered 
state of our system cannot be globally stable, e.g. [56, sec. III]. Since the energy 
develops huge values (£7 « 10 6 at p B = // = 2M q fig. |24"|) . this also suggest that 
the considered state is globally unstable. 

In QCD, it is expected that mesons condense at small temperature and large 
isospin chemical potential, which corresponds to our unstable region. See [57-59] 
for examples of vector condensation and [46,60,61], as well as [62,63] for examples 
of pion condensation. We expect similarly to QCD a new stable phase, which is 
determined by the condensation of mesons. These mesons are dual to fluctuations 
about the background. In section [6] we will discuss these fluctuations and indeed 
find the expected new phase. 



5 Hydrodynamics 

In this section we consider the effective diffusion coefficient computed from the 
membrane paradigm developed in [10] and extended in [39]. The diffusion coeffi- 
cient for the pure baryonic case is presented in [43] . To include isospin density we 
use the following procedure. The finite baryon and isospin density enter the dif- 
fusion coefficient through the D7 embedding function x(p>" ,cr) which appears 
in the metric G(p,d B ^d 1 ) A2.6H . We obtain the explicit embedding function by 
solving its equation of motion (|2.43p and then simply plug in the metric into the 
diffusion formula for the case of zero densities^ 



D 



G11V— G00G44 



P =i 




(5.1) 



The calculation of D makes use of the isospin depending embedding function 
x{p,d B ,d I ). Thus, we call it the effective baryon diffusion coefficient. This pro- 
cedure yields the plots given in figure [25] (a) and (b) . 

In the limit m — > 00 and m — > 0, the diffusion coefficient approaches always 
D = l/(27rT), which coincides with the result for the diffusion coefficient of R- 
charges in N = 4 SYM. For zero densities the diffusion coefficient D vanishes 
at the phase transition because all the quarks are bound to mesons [39]. If we 
increase the baryon and isospin densities, there are always free quarks available and 
the diffusion coefficient never vanishes [43]. For densities larger than the critical 
one, the first order phase transition disappears. The only structure in DT which 
survives in the large densities region is a minimum. We consider this minmum 
of the diffusion coefficient as a hydrodynamic crossover point. Including these 
crossover points in the thermodynamical phase diagram from fig. 1141 (a), we obtain 



13 Notice that a strict dervivation of the diffusion coefficent for the case of finite densities would 
change equation (15. ip such that terms explicitly depending on the finite densities appear. We 
leave this task for further studies and restrict our analysis to equation Q5.10 . 
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Figure 25: The diffusion coefficient DT versus the mass parameter m as defined in 
equation (|2.45p at d B — 5 • f -5 (a) and d B = 20 (b). The five different curves correspond 
to d 1 = (black), d 1 = \d B (green), d B = \d B (blue), d B = \d B (red) and d 1 = d 1 
(orange) . (c) is the " enlarged" phase diagram, where we plot the position of the minimum 
of the diffusion coefficient versus the densities, (d) shows a contour plot of the "enlarged" 
phase diagram, which demonstrates the breaking of the rotational symmetry 0(2) down 
to the discrete symmetry Z4. 
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fig- El] (c), where the surface marks the phase transition and cross over point, 
respectively. As we expect, the rotational symmetry 0(2) is completely broken 
for large densities and in the case Nf = 2 there is only the discrete symmetry Z4 
left. 

It is important to notice is that the baryon density increases the diffusion 
coefficent as we see in fig. [25] (a), (b) and [43]. Increasing the isospin density 
instead reduces the diffusuion coefficient (see fig. [25] (b)). Therefore, the effects 
of the baryon and the isospin density on this effective baryon diffusion coefficient 
are opposite. 

6 Mesons at finite chemical potential 

In this section we use the AdS/CFT correspondence to calculate the mesonic 
spectral function of the field theory by investigating the fluctuations of D7-branes. 
We then study the effects of chemical potentials on the meson spectrum. 

As in [39] , we interpret resonances in the spectral functions of flavor currents as 
quasiparticle bound states, which are identified with mesons. A direct verification 
of this identification was given in [43]. There the congruence of the spectra ob- 
tained from spectral functions and the analytically known meson spectra from [44] 
was shown in an appropriate limit. 

6.1 Mesonic bound states and brane fluctuations 

To describe mesons, we compute the two point correlation function of flavor cur- 
rents J(x) of the field theory. The retarded Green function G R in momentum 
space may be written as 

G R {u,q) = -i yd 4 xe ife 0(x°)([J(x),J(O)]). (6.1) 

We write k for the four-momentum vector with spatial three-momentum q and 
timelike component 00, such that k = (u, q) T . The description of the spectrum of 
mesons is given in terms of the spectral function d\(uj, q), defined as 

m(cu,q) = -2lmG R (u,q). (6.2) 

Poles in the Green function give rise to peaks in the spectral function, which we will 
study below. The position, magnitude and width of these peaks are determined by 
the position and structure of the poles of G R . They encode the mass and lifetime 
of the mesons [36-38] . 

To compute the current correlator in (|6.1|) , we use the AdS propagator of the 
relevant supergravity field. The dual gravity field to the scalar, pseudo-scalar and 
vector flavor currents are the fluctuations of the probe D7-branes and the fluctu- 
ations of the gauge field on the branes. In terms of the coordinates introduced 
in section [2] the scalar currents are dual to fluctuations of the scalar embedding 
function 9(p), the pseudo-scalar currents correspond to fluctuations of the angular 
embedding (j)(p), and the vector current is mapped to fluctuations of the gauge 
field A{p) on the branes. 
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6.1.1 Calculation of spectral functions 



The computation is performed following the prescription of [34,35]. Moreover, 
in this section we extend the results of [43] and make use of the calculations 
performed there. 

The method of calculating spectral functions in a holographic model amounts 
to solving the equations of motion for a supergravity field A(p, q) in AdS space. 
Then evaluate 

G H ( „, g) = w! lim 



p—rOO 



A(p,k) 



(6.3) 



to obtain the Green function and eventually make use of (16. 2p . which gives the 
spectral function in terms of the solution A(p, k), 

9^,S)—M^. jBl t nl U?tM). (6.4) 
4 p^oo I A(p) J 

Note that the formulae for G R and d\ above describe the result of the procedure 
described in [34], tailored to the field and brane configuration used in this paper. 

The remaining task is to derive and solve the equations of motion for the field 
A. We will now follow this procedure for the cases of non-vanishing baryonic 
or isospin chemical potential in the case of a theory with Nj = 2 flavors. In 
particular we concentrate on the vector mesons, dual to fluctuations A of the 
gauge field A on the D7-brane. For simplicity, we restrict ourselves to the limit 
of vanishing spatial momentum q = 0. The vector mesons show the interesting 
feature that the spectral lines split up in the case of non-vanishing isospin chemical 
potential [43]. In section we will make use of this property to explore the phase 
diagram of fundamental matter in the D3/D7 setup. 



6.2 Spectra at finite baryon density 
6.2.1 Equations of motion 

In the case of pure baryonic chemical potential the relevant action is the DBI- 
action (I2.16p . which now takes the form 

S = -T D7 J d 8 £ Str vldet (G + lixa' F)\. (6.5) 

We have introduced chemical potentials in section I2T21 and derived the equation of 
motion for the background gauge field A and their solution in section 12.3.21 

To calculate the mesonic spectral functions, we now investigate small fluctu- 
ations of the background field configuration, which extremizes the action (|6.5|) . 
Therefore we add small fluctuations A to the background A, 

A^A + A, (6.6) 
F^F + F. (6.7) 

Note that according to section 12.21 only the timelike component A${p) of the 
background field is non-vanishing. The fluctuations, however, are gauged to 



45 



have non-vanishing components along all Minkowski directions and depend on the 
Minkowski coordinates x and on the radial coordinate p, -Aoi,2 3 = ^0,1,2,3(^/5) 
and A4, = 0. We write 

g = G + 2na'F (6.8) 
and in this way cast the Lagrangian for the fluctuations into the form 

(6.9) 



det [g + 2vra' F 



C = 

In general the non-linear structure of the Lagrangian induces couplings between 
vector mesons, scalar and pseudo-scalar mesons, as there will be couplings of 
fluctuations in Q and those of F = dA. However, in the limit of q — > 0, which we 
are working in, these couplings do not occur [42]. We assume the fluctuations A 
to be small, justifying the consideration of the linearized equations of motion for 
the field components A^, 

v^detgj (g^g a '< - g^g u '< - gWff*A aij , (e.io) 



o = d v 

where we write upper indices on g to denote elements of Q . 

As in [39] we should consider gauge invariant combinations of field compo- 
nents only. For arbitrary momentum we are free to choose a reference frame in 
Minkowski space, such that k = (cj,q, 0, 0), and we have gauge invariant longitu- 
dinal and transversal components 

E\\ = loA x + qA , E ± = uA x>z . (6.11) 

In the limit q — > 0, the equations of motion (and therefore also the according 
Green functions) for all components become identical, since there is no distinction 
between longitudinal and transversal modes. 
In momentum space, 

MP, x) = j-^ e ifc %(p, k) , (6.12) 
the equations of motion (|6.1U|) for the gauge invariant fields are given by 



^det^G 22 G 44 G u 
- E" + i),, Id ( P 3 f(l~X 2 ) 2 I L < (G.13) 

/l_ x 2 +p 2 x ,2_ 2 j_ {1 _ x 2 ){dpAo) 2^ 

with dimensionless It) = ui/(2irT), and dimensionless coordinate p as well as di- 
mensionless fields A \—* ^^-A. A prime denotes a derivative with respect to p. 
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Figure 26: A spectral function oscillating around the asymptotic result EHo- We plot 
spectral functions in units of N fN c T 2 / '4. 



6.2.2 Numerical results and interpretation 

We solve the equation of motion (16.130 numerically, using in-falling wave boundary 
conditions at the black hole horizon. Subsequently we make use of (j6.4j) (where 
A i— ► E) to compute the spectral function. Remember that this function depends 
parametrically on the particle density (or on the chemical potential, respectively), 
which influences the background field solutions A and the embeddings x in fl6-13|) • 
The large to behavior of the spectral functions can be found analytically and 
is given by 

9to(tt),0) = N f N c T 2 7TXo 2 . (6.14) 

An example for a spectral function and the oscillations around 9\q can be found in 
fig. [26] One clearly sees the various excitations as peaks in the spectral function. 
Each peak corresponds to a quark-antiquark bound state with vanishing angular 
momentum. 

The position ro n of a peak corresponds to a (dimensionless) meson mass M n = 
\k\/(2nT) = tt) n . In [43] it was observed that the meson masses obtained in this 
way indeed agree with the spectrum from [44], found analytically in the super- 
symmetric zero temperature limit d — > 0, m —* oo, 



M n / (w + 2)(n+l) 

2^T =m V 2 ' ri = 0,1,2,..., (6.15) 

where m is the parameter defined in (j2.45jl . At finite temperature and baryon 
density we will observe deviations from this spectrum. Nevertheless, at small 
densities and not too small m, the first peaks in the spectrum still come very close 
to this formula, see fig. [26] with the first peak corresponding to n = 0, the second 
to n = 1 and so on. 

The dependence on temperature is related to the dependence on the quark 
mass, since the relevant parameter in the D3/D7-setup is the quotient m. The 
dependence on this parameter was intensively investigated in [43], extensions to 
finite momentum and the coupling to scalar and pseudo-scalar modes can be found 
in [41,42]. 
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Figure 27: Spectral functions for various baryon densities d, again normalized to 
NfN c T 2 /4:. The black and brown curves for small d show meson states around roo = 5 
and higher excitations. At higher values of d these excitations disappear, and at very high 
densities a new structure forms at small ro. 

We work in the canonical ensemble and will now investigate the effects of 
variations in d. Spectral functions for various finite baryonic d are shown in 
fig. [53 Again, the peaks in these spectral functions indicate that quarks form 
bound states. At low baryon densities the positions of the peaks agree with the 
supersymmetric result (|6.15p . e. g. the peaks at roo ~ 5 and roi as 8.7 in the black 
curve correspond to the excitations with n = and n = 1 in (I6.15|) . Increasing 
the quark density leads to a broadening of the peaks, which indicates decreasing 
stability of mesons at increasing baryon density [48,64]. At the same time the 
positions of the peaks change, which indicates a dependence of the meson mass on 
the baryon density. Further increasing the quark densities leads to the formation 
of a new structure at to < 1. We will discuss this structure together with the 
results at finite isospin density. 

6.3 Spectra at finite isospin density 
6.3.1 Equations of motion 

To investigate isospin effects, we again start from the action (|2.16p . This time we 
choose F as in ([2TTUT) with only F 1 ^ 0. Since the /9-dependent time component of 
the gauge field generates the chemical potential on the field theory side, again we 
only have non-zero components Fg = —Fq 4 in the background. Now a is an index 
labeling the SU(2) components. So our gauge fields A not only carry Lorentz- 
indices but additionally SU(2) gauge indices, as in A a ^ with a = 1, 2, 3. In general 
the field strength tensor therefore has elements F® v = 2d^A^ + f abc A b ^A^, with 

structure constants f abc = e abc . 

In the non-Abelian field strength tensor, the term quadratic in the gauge field 
describes a self interaction of the gauge field. The coupling constant for this 
interaction may be determined by a redefinition of the gauge field, such that the 
kinetic term of the effective four-dimensional theory has the canonical form. In 
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appendix [D] we show that the redefinition is given by 



A^^A, (6.16) 

where the dimensionless constant c depends on the geometry of the D7 worldvol- 
ume directions along p and on the S 3 , which are transverse to the directions of the 
D3. In particular, c is independent of the 't Hooft coupling A. Determining the ex- 
act value of c is left to further work in terms of the ideas presented in appendix [DJ 
In the following we set c = 77=. The field strength tensor in the redefined fields is 
given by 

F« v = 2d^A a v] + ^f ahc A\Al (6.17) 

We make use of the derivations performed in [43] (where the above redefi- 
nition was implicitly used) and assume the isospin chemical potential to be ori- 
ented the direction labeled by a = 3, as motivated in (12.100 . Therefore the only 
non- vanishing background component is -F 3 = — i*g 4 , which coincides with the 
background for the case of a baryonic chemical potential. 

In dimensionless quantities, as in the previous subsection, the linearized equa- 
tions of motion in the case of vanishing spatial momentum for the components A't 
of gauge field fluctuations A can be derived as 



= d K 



V\detg\ {Q VK Q°» - G VU G^) F 



(6.18) 



2ira 



where P", = 2d^A$, + c/\/A/ a63 ^(<5o M ^ + Q 2 H /{2ira') contains all field 

strength contributions linear in A. 

To solve these equations of motion, we first introduce the gauge invariant fields 
E a as in (16. lip and then perform another change of basis to express the equations 
of motion in terms of the fields 

X = E 1 + iE 2 , Y = E 1 -iE 2 , E 3 . (6.19) 

In this basis the equations of motion decouple and may be written as 



= X " + PVV l-x' -A^—r (tv-m) 2 X, (6.20) 

^detgie? 44 ^ 22 r a q m v ' y ' 



d p MdetS e 44 S 2i ,4 goo 

= Y" + V , >-Y> -4%^ (w + m) 2 Y, (6.21) 

^detgje? 44 g 22 R A Q iA 



= E 3 + — \ J -E 3 -4^r^—to 2 E 3 , 6.22 

y|detg|g 44 g 22 R A Q M v 

where we introduced 
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Figure 28: Spectral function calculated from (16.211) . A new peak forms at high densities. 



Here -Aq(p) ^ s t ne dimensionless background gauge field. The quantity m is related 
to the dimensionful background field Aq (p) by m = • 

As a result of the introduction of an isospin chemical potential instead of a 
baryonic one, we end up with three equations of motion for the decoupled fields X, 
Y and E 3 . These equations now will have three distinct solutions because of the 
differences in the last terms of (16.201) to (I6.22f) . Since the difference is a shift in ro 
one would expect the solutions X and Y to be shifted versions of E 3 . However, the 
shift parameter m does depend on p. The equation of motion for E 3 is identical 
to the baryonic case and as such was discussed in the previous subsection. 

The three solutions X, Y and E 3 constitute the isospin triplet of mesons 
which may be constructed out of the isospin 1/2 quarks of the field theory. This is 
similar to the p- meson in QCD. The mode E 3 coincides with the solution in case 
of a pure baryonic chemical potential, while the other two solutions have peaks in 
the spectral function at lower and higher values of ro [43]. The magnitude of this 
splitting of the spectral lines is determined by the chemical potential. 

In the limit of zero frequency ro — * 0, equations (|6,2U|) and (|6.21|) coincide and 
will result in identical solutions X and Y. In this limit the solution E 3 , though, 
differs from X and Y, by means of the last term. So for small frequencies to, we 
expect differences between the solutions E 3 and X, Y. 



6.3.2 Numerical results and interpretation 

We now investigate the effects of finite isospin density on the spectrum. After 
solving (|6.20p to (|6.22p numerically, the correlators may eventually be evaluated 
as given by (|6.4p . where A has to be replaced by E 3 , X or Y. The peaks in the 
spectral functions again correspond to mesons. 

An interesting feature at finite isospin chemical potential is the formation of a 
new peak in the spectral function in the regime of small ro at high density/high 
chemical potential, see fig. [5HJ Notice that compared to the baryonic case the 
density, at which the new peak forms, is about two orders of magnitude smaller. 
As in the baryonic case, the excitations related to the supersymmetric spectrum 
broaden, the corrsponding mesons become unstable. 
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Figure 29: Plot of the spectral function around in = 0. At a value of d = 15.35 a pole 
appears at ro = 0. 

One should keep in mind that the structure of the spectral function is deter- 
mined by the pole structure of the retarded correlator (|6.1[) . The poles of this 
function are located in the complex w-plane at positions f2 n £ C. The spectral 
functions show the imaginary part of the correlator at real valued uj. Any pole 
in the vicinity of the real axis will therefore introduce narrow high peaks in the 
spectral function, while poles far from the real axis have less influence and merely 
introduce small and broad structures. 

In reminiscence of the real valued normal modes of solutions for oscillators the 
complex Q n are called quasinormal modes. The imaginary part of the quasinormal 
modes describes damping, as long as Imf2 n < 0. The figures shown in this section 
indicate a dependence of the position of the quasinormal modes on the chemical 
potential/particle density. From fig. [28] we deduce that a quasinormal mode ap- 
proaches the origin of the complex u plane as the particle density is increased. 
We observe a pole at to = for a certain particle density do-it, the value depends 
on m. An impression of the variation in the spectral function is given in fig. [53 

In fig. [30] we qualitatively sketch the result from the investigation of the be- 
havior of the quasinormal modes closest to the origin of the complex ro-plane. 
These modes do not produce the peaks corresponding to the spectrum (|6.15p . At 
low densities all quasinormal modes are located in the lower half plane. When 
increasing the isospin density, the modes of the solutions X and Y to f|6.20|> and 
(|6.2ip move towards the origin of the frequency plane. At the same time two 
quasinormal modes of E 3 move towards each other and merge on the negative 
imaginary axis, then travel along the axis towards the origin as one single pole. 
At the critical value of d the modes from X and Y meet at the origin, the quasinor- 
mal modes from E 3 still reside in the lower half plane. This observation matches 
the discussion below equations (|6.2U|) to (|6.22|) . where we expected X and Y to 
behave similarly at small to, while E 3 should differ from this behavior. Upon 
further increasing the isospin density, the modes O from X and Y enter the upper 
half plane, maintaining their distinct directions. The sign change in Imf! from 
Imfl < to ImJ7 > indicates that a damped resonance changes into a self- 
enhancing one, and thus introduces an instability to the system. Figure fig. [31] 
illustrates the transition of a quasinormal mode of Y from the lower half plane to 



: d = 14.6 

- d= 15.8 




m — 3 
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Figure 30: A sketch of the positions and movements of the quasinormal modes under 
changes of d. Color indicates the function: red = Y, green = A, blue = E 3 . The symbols 
indicate the range of d: o < c? cr jt, • = d cr n, ■ > ci C rit- Poles in the gray region introduce 
instabilities. 
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Figure 31: Contour plot of the spectral function for Y around tt> = in the complex 
tt>-plane. Left: d = 10 < d C rit, center: d = 15.352 = d C rit 5 right: d — 20.704 > d CY i t - The 
three graphs where generated for m — 3, dark shading indicates small values of £H, light 
shading indicates large values. A pole in the upper half plane introduces an instability. 



the upper half plane. The E 3 -mode does not enter the upper half plane at any 
value of d we considered. Compare this to the values of d in fig. [27] and fig. [28] 
at which the pole induces visible structures at small ro. A comparable movement 
of poles in a different but related setup was found in [65]. There the quasinormal 
modes of correlation functions of electromagnetic currents were investigated as a 
function of temperature. 

In the following we interpret the observation of decaying mesons and the emer- 
gence of a new peak in the spectral function in terms of field theory quantities. 
In particular we speculate on a new phase in the phase diagram for fundamental 
matter in the D3/D7 setup. 

In the far UV, the field theory dual to our setup is supersymmetric, thus 
containing scalars as well as fermions, both of which contribute to the bound 
states we identified with mesons, even when supersymmetry is eventually broken. 
The meson decay at non-vanishing particle densities may be explained by the 
change of the shape of the potential for the scalars in the field theory upon the 
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introduction of a non-vanishing density. As outlined in appendix [Cl a chemical 
potential may lead to an instability of the theory, since it induces a runaway 
potential for the scalar fields at small field values [66]. Nevertheless, interactions 
of 4 -type lead to a Mexican hat style potential for larger field values. In this way 
the theory is stabilized at finite density d while the scalar fields condensate. This 
squark condensate presumably contributes to the vev of the scalar flavor current, 

d oc <J°> oc ($7°V) + (4>d G 4>) . (6.24) 

In the AdS/CFT context the presence of an upside-down potential for the squark 
vev has been shown in [50] using an instanton configuration in the dual supergrav- 
ity background. 

The occurrence of a pole in the upper half plane of complex frequencies at 
finite d CT - lt indicates an instability of the theory. A comparable observation was 
made in [48], where in fact the vector meson becomes unstable by means of neg- 
ative values for its mass beyond some critical chemical potential. The difference 
between this work and [48] is that our model includes scalar modes in addition 
to the fundamental fermions. Nevertheless, in both models an instability occurs 
at a critical value of the chemical potential. The theory may still be stabilized 
dynamically by vector condensation [67]. In this case the system would enter a 
new phase of condensed vectors at densities larger than d CT \t, in accordance with 
the expectation from QCD calculations [57-59]. 

We perform the analysis of the pole structure at in = for various m, and 
interpret the phenomenon of the transition of poles into the upper half plane at 
finite critical particle density as a sign of the transition to an unstable phase. We 
relate the critical particle density d cr - lt to the according chemical potential p, 1 by 
(|2.37p and use the pairs of m and critical dimensionful /i to trace the line of the 
phase transition in the phase diagram of fundamental matter in the D3/D7 setup. 
The result is drawn as a green line in fig. [32l The picture shows the (/x 7 , T)-plane 
of the phase diagram and contains three regions, drawn as blue shaded, white, and 
green shaded, as well as a blue and a green line, separating the different regions. 

The blue shaded region marks the range of parameters, in which fundamental 
matter is described by D7-branes with Minkowski embeddings. The blue line, 
delimiting the the blue region, marks the line of phase transitions to the black 
hole phase, where fundamental matter is described by D7-branes which have black 
hole embeddings. This phase transition was discussed in section 14.21 Using the 
symmetries ()2.33j) . this phase transition line can be mapped to the line of phase 
transitions between Minkowski and black hole embeddings, present at finite baryon 
chemical potential [28,51,68]. 

The green line and green shaded region in the phase diagram in fig. [32] mark 
the observation made in this section. The green line marks the values of d cr a at 
which the pole in the spectral function appears at tt) = 0. Beyond the green line 
we enter the green shaded unstable region. 

We observe that the green line asymptotes to a straight line at high temper- 
atures. Within the values computed by us, this line agrees with the asymptotic 
behavior of the contour of particle density with d ~ 20.5, drawn as a gray line 
in the phase diagram. We thus speculate on a finite critical particle density be- 
yond which the black hole phase is unstable. This interpretation is supported by 



53 





3.0 r 




2.5 1 




z.U - 








i.5 ; 


a. 






1.0 




o.5 ; 




0.0 1 











d = 20.5 




unstable 








black hole 




Minkowski 







0.0 



0.2 



0.4 0.6 

1/m oc T/M q 



0.8 



1.0 



Figure 32: The phase diagram for fundamental matter in the D3/D7 setup. In the blue 
shaded region D7-branes have the topology of Minkowski embeddings. The white and 
green regions are described by black hole embeddings. However, in the green region we 
interpret our observations as an instability of the black hole embeddings in the D3/D7 
setup. 



analogous studies of the phase diagram of N = 4 super- Yang- Mills theory with 
R-symmetry chemical potentials, where a similar line in the phase diagram was 
discovered [45,56]. Up to now it is unclear whether the new phase indeed is unsta- 
ble or if it only indicates that black hole embeddings are not capable of describing 
this state of matter. 

Note that the location of the green phase transition line in fig. [32] as well as 
the results shown in fig. [31] and fig. [30] are obtained from the analysis of poles 
in the spectral functions. These functions in turn are obtained as solutions to 
equations (|rj.2U|) to (|6.22|) . which do depend on the so far unknown factor c in 
(|6.23p . The computation of this factor is left to future work. It will determine the 
exact position of the green phase transition line in fig. [3"2"1 This will answer the 
question whether there is a triple point in the phase diagram and if the blue and 
green shaded regions meet at a common border. Moreover, other poles than the 
ones investigated here may have influence on the stability of this system. 

7 Conclusion 

In this work we have considered a probe of Nf D7 branes embedded in the AdS 
Schwarzschild black hole background with an U(Nf) gauge field. This configura- 
tion is conjectured to be dual to a thermal quantum field theory with isospin and 
baryon charge. We have considered both the canonical and grand canonical en- 
semble. Our detailed study of this background provides a basis for further studies. 
For instance it is now possible to investigate the mesonic excitation spectrum or 
hydrodynamic properties in detail. 

By including isospin density, we find four new phase boundaries coincident 
with the fixed points of the accidental symmetries between baryon and isospin 
density, which our setup displays. Along the lines of equal densities \d B \ = \d J \ 
and near T = 0, we find a sharp transition in the chemical potentials. This 
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situation is reminiscent of a corresponding transition found in two-color QCD [46]. 
For T = the chemical potential is discontinuous and a one-to-one correspondence 
between the thermodynamic conjugate variables d and p, can not be established. 
This inconsistency between the two ensembles signals an instability of the setup 
in that phase region. We expect that vector mesons condense in this region as is 
also expected in QCD [46]. A further detailed study of this new phase requires 
the investigation of scalar and pseudoscalar fluctuations as well. Furthermore, an 
analysis of the possible condensates is required, cf. [48] . This is a promising task 
for future investigation. It would also be interesting to study the heat capacity 
and the quasi-particle spectrum at low temperature as in [69] in the more general 
setup provided in this paper. This study may give a better understanding of the 
new quantum liquid found in [69] . 

Here we have restricted to vector fluctuations about the finite isospin density 
background. Analyzing these fluctuations, we found an unstable region in the 
(/u 7 , T) phase diagram which is bounded by a line of constant density. We expect 
that mesons condense in this unstable phase. The emerging phase structure ap- 
pears to be a generic feature of holographic duals with chemical potentials as a 
similar behavior has been observed for spinning branes [70] and R-charged black 
hole backgrounds [56,71]. It would be interesting to shed more light on the connec- 
tion between the new phase found in this work and in the backgrounds mentioned 
above. 

Considering the effective baryon diffusion, we find competing behavior between 
baryon and isospin density. For example at a given temperature, increasing d B in- 
creases the effective baryon diffusion while increasing d 1 suppresses that diffusion. 
Note that this effective baryon density collects both, the effect of finite baryon 
and of finite isospin density on the effective baryon diffusion. It would be very 
interesting to separate these two contributions, for example by a modification of 
the membrane paradigm. 

All the numerical results mentioned here were obtained for Nt = 2. However, 
they may also be generated for arbitrary Nf <C N c from the general expressions 
we have provided here. 
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A Simplification of the DBI action 

In this section we give some details for the simplifications done in the equations 
(|ZZE)) , flZ22D , ([223]). We argued that the square root in the action (|2T2Tj) is 
diagonal in flavor space and the DBI action decouple for the fields defined in 
(|2.22p . To prove this, we first have to calculate the square of the non- vanishing 
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field strength component (F^q) 2 . For this calculation we need some multiplication 
relations for the matrices A* given by 



,...,N f 




(A 1 ) = diag(l, . . . , (JV> - 1)', . . . , 1) i = 2,...,N f 



(A.l) 



i-th position 

A^' = diag(l,.../-(iV/-l),. 



j'-th position 

H/V/^T), . . . , 1) i + j 



and the fact that all matrices X 1 commute with each other since they are the 
generators of the Cartan subalgebra. The i-th diagonal component of the matrix 
(-F40) 2 is then for i = 1 



((^40) 2 ) 



11 



N f N f N f N f 



3=1 



3=2 



jr'=2 fc=2 



J 11 



v3=l 



10 



(^) 2 



(A.2) 



and for « 7^ 1 



3=1 
3 A 



£(i^) 2 (A^') 2 + (F 4 ^) 2 (V) 2 + 2^F^AV 

3=1 
3 A 

+ 2 EE i? 4>4 / o fc A J A fc 
/ 



j=i fc=i 



(A.3) 



\ 



$>2 - (JV> - 



\ J=1 



J 



where we used the new fields as defined in equation (12.220 in the last equalities. 
With this choice of the new fields and since -F40 is diagonal in the flavor space, 
each component of the diagonal is given by the expression known from the pure 
baryonic case [27] 



A 1 + (2vra') 2 G' 00 G 44 ( J F4o) 2 .. = \A + (27ra / ) 2 G 00 G 44 ((F 4 o) 2 ) i 



l + (2vra / ) 2 G 00 G 44 (^X i ) 2 . 



(A.4) 



Taking the trace and simplifying the expressions of the metric factors, we get 
equation (|2 .231) . 
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B Effect of the accidental symmetries for Nf = 3 



In thi s section we give the action of the accidental symmetries discussed in section 
12.3.11 in the physical basis for Nf = 3: 





d 2 : 


d B 


^d B 




d'2 


I-+ -(d'2 + d J 3) 


d'3 


^d J 3 


d! « 


d 3 ■■ 


d B 


~d B 




d'2 


^d J 2 


d'3 


I-+ -(d 7 2 + d J 3) 


d 2 « 


d 3 : 


d B 


~d B 




d 1 ? 


^d J 3 


d'3 


^d J 2 


d! « 


-di : 


d B 


1 H 
^-,d B 


- 2d J 2 _ 2d 7 3 ) 


d 1 ? 


w ^(- 2dB + d' 2 - 2[jl3 ) 


d T 3 


w -(-2d B - 2d'2 + d 7 3 
3 


d 2 « 


-d 2 : 


d B 


^-,d B 


+ 2d 7 2) 


d 1 ? 


I — * ^(4d B - d 7 2) 


d'3 


w -(-2d B + 2d J 2 + d 7 3 
3 


d 3 


-d 3 : 


d B 




+ 2d'3) 


d 1 * 


_ ^(-2d B + d' 2 + 2d 7 3) 


d J 3 


_ i ( 4d B -d 7 3). 

3 



Since the transformation matrix to the physical basis for Nf = 3 is not an or- 
thogonal matrix, there is no induced approximate 0(3) symmetry in the physical 
basis. This is in general also true for all cases with Nf > 3. 



C Chemical potentials in field theories: Runaway po- 
tential and Bose-Einstein Condensation 

In our setup we consider a theory which is super symmetric in the far UV. Its 
fundamental matter consists of complex scalars (squarks) and fermionic fields 
(quarks). In this section we describe the effect of the chemical potential on the 
field theory Lagrangian and on the vacuum as e. g. in [66] . We consider a theory 
with one complex scalar (j) and one fermionic field tp with the same mass M q 
coupled to an U(l) gauge field A u . The time component of the U(l) gauge field 
has a non-zero vev which induces the chemical potential fi, 

A u = fi5 u0 . (C.l) 

The Lagrangian is given by 

C = -D v <pD v 4> - M 2 q 4>*4> - i;(D+M q )4, - , (C.2) 

where D u = d u — \A V is the covariant derivative and F^ u = d^A u — d u A^ the field 
strength tensor. Expanding the Lagrangian around the non-zero vev of the gauge 
field, the Lagrangian becomes 

£ = -d u p*d u <p - (M 2 - ft 2 )<j)*<j) + ^J S - + M q )^ + /iJ F , (C.3) 

where = i[(c?^*)0 — <£*(c^0)], (J^) F = —iip^^ip are conserved currents. These 
conserved currents are the population density for the scalar N$ and fermionic 
fields Nf, such that the linear terms in the Lagrangian are /jlNs and fiNp- 

The mass term — {M 2 — ^ 2 )4> 2 of the Lagrangian (|C.3|) introduces an insta- 
bility if fj, > M q since the corresponding potential V = (M 2 — fj, 2 )(j) 2 + • • • is 
not bounded from below. In some systems this runaway potential is stabilized 
by higer interactions and becomes a Mexican hat potential such that the scalar 
condense and the scalar density becomes non-zero. This condensation is known 
as Bose-Einstein-Condensation. 
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D Coupling constant for vector meson interaction 



In this section we show how the coupling constant for the interaction of vector 
mesons can be computed, extending the ideas presented in [44]. This coupling 
constant in the effective four-dimensional meson theory can be determined by 
redefinition of the gauge fields such that the kinetic term has canonical form. 
This coupling constant depends on the geometry of the extra dimensions. 

First we consider the eight-dimensional theory determined by the DBI action 

(2) 

£>dbi expanded to second order in the fluctuations A, 

41 = TD7(2 4 TO/)2 / G^G^F^F^ , (D.l) 

where G contains the background fields and we simplify the analysis by considering 
only Abelian gauge fields. Defining the dimensionless coordinate p = g/R and 
integrating out the contribution of the £ 3 , we obtain 

^ = ^(WNW I d 4 x J dpy/ ^ g^g^F^F^ . (D.2) 

To obtain a four-dimensional effective theory we have to integrate over the coordi- 
nate p. This contribution depends on the geometry induced by the p dependence 
of the metric factors. However, we expect that it is independent of the 't Hooft 
coupling A. We parametrize this contribution by d . The kinetic term of the 
effective theory is then given by 



where the pref actor may be written as 

T D7 (27ra / ) 2 vol(5 3 )« 4 c / 



> (D-4) 

9ym c 



where the numerical values independent of the 't Hooft coupling are grouped into 
the coefficient c. From this we can read off that a rescaling of the form 

A h-> -^=A (D.5) 
casts the Lagrangian into canonical form with a prefactor of 1 /<7ym • 
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